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INTRODUCTION 


This  report  covers  progress  made  during  the  past  year  in  the  work  of  eleven  Re¬ 
search  Units  funded  under  the  Joint  Services  Electronics  Program  at  Harvard  University. 
It  is  broken  down  into  four  major  divisions  of  electronic  research:  Solid  state  electronics , 
Quantum  electronics ,  Information  electronics ,  and  Electromagnetic  phenomena.  Following 
the  report  of  the  work  of  each  Unit,  there  is  a  complete  annual  report  of  the  associated 
Publications/Patents/Presentations/Honors.  This  report  also  includes  a  section  on  Signif¬ 
icant  Accomplishments  which  contains  selected  highlights  from  two  of  these  Units  entitled 
“Competitive  Analog  Neural  Networks  for  Translation-Invariant  Feature  Detection,”  Re¬ 
search  Unit  5,  and  “A  Novel  Resonant  Array  of  Parallel  Dipoles  with  Unique  and  Useful 
Properties,”  by  Research  Unit  11. 
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1.1  Electronic  Theory  of  Semiconductor  Alloys  and  Superlattices.  H. 

Ehrenreich,  Grant  N00014-89-J-1023;  Research  Unit  1. 

During  the  past  year  the  theoretical  research  efforts  of  this  group,  consisting  of 
C.  Grein,  P.M.  Hui,  E.  Runge,  P.M.  Young  and  H.  Ehrenreich,  have  focussed  on  1)  ex- 
citons  and  interband  transitions  in  III-V  semiconductor  superlattices,  2)  the  response 
and  transit  times  in  quantum  well  structures,  3)  phonon  assisted  transport  through 
double  barrier  resonant  tunneling  structures  and  4)  Auger  processes  in  GalnSb/InAs 
and  HgCdTe  IR  detectors.  This  year’s  research  output  in  terms  of  published  papers 
was  diminished  because  of  H.  Ehrenreich’s  involvement  with  the  new  Federal  Program 
on  Materials.  Ehrenreich  spent  much  of  the  fall  1991  at  the  White  House  as  Consul¬ 
tant  to  Presidential  Science  Advisor  D.A.  Bromley  and  helped  to  formulate  plans  for 
implementing  this  initiative.  At  Bromley’s  request  he  prepared  a  white  paper  entitled 
The  Proposed  Federal  Program  in  Materials  Seienee  and  Technology:  Comments  and 
Views  (1991,  unpublished). 
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A.  Excitons  and  Interband  Transition*  in  ill-V  and  II- VI  Semiconductor  Superlat¬ 
tices 

The  fundamental  absorption  in  InGaAs/InAlAs  and  GaAs/GaAlAs  superlattices 
has  been  calculated  quantitatively  using  a  superlattice  K  •  p  theory.  Electron-hole 
Coulomb  interactions  yielding  excitons  and  interband  transition  enhancement  associ¬ 
ated  with  the  Sommerfeld  factor  are  incorporated  in  the  theory.  The  Coulomb  inter¬ 
action  is  treated  within  an  independent  subband  approximation  that  associates  each 
exciton  with  a  single  conduction  and  a  single  valence  band.  This  approach,  together 
with  a  simple  model  representing  the  superlattice  Wannier  functions  which  is  both 
simple  and  accurate,  leads  to  an  analytic  form  for  the  electron-hole  interaction  that 
permits  calculation  of  the  exciton  wavefunction,  binding  energy  and  optical  absorp¬ 
tion  without  resorting  to  a  variational  approach.  The  only  input  consists  of  constituent 
bulk  parameters,  band  offsets  and  strains  as  described  by  deformation  potentials.  The 
calculations  require  only  10  minutes  of  CPU  time  on  a  SUN  system,  and  yield  optical 
structure  within  2-3  meV  and  absolute  absorption  coefficients  within  10%  of  experi¬ 
mental  results  for  all  but  one  of  the  12  InGaAs/InAlAs  and  GaAs/GaAlAs  samples 
analyzed.  The  experimental  optical  absorption  data  and  the  sample  quality  chosen 
for  comparison  with  the  theory  are  believed  to  be  the  most  reliable  available  at  the 
present  time.  The  modest  computer  requirements  permit  an  investigation  of  the  behav¬ 
ior  of  optical  structure  as  a  function  of  band  offsets  and  other  interface  parameter  and 
therefore  allows  estimates  of  physically  important  quantities  of  relatively  unexplored 
heterostructure  systems,  for  example  II- VI  superlattices  [1]. 

This  same  approach  has  been  applied  to  wide  band  gap  II- VI  superlattices  of 
possible  importance  to  blue-green  laser  applications.  The  excitonic  and  interband 
optical  absorption  in  ZnCdSe/ZnSe  and  ZnCdTe/ZnTe  have  been  explored  utilizing 
the  data  of  A.  Nurmikko’s  group  at  Brown  University.  A  paper  in  preparation  describes 
the  results  of  this  work  [2].  The  II- VI  materials  are  distinguished  from  their  III-V 
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counterparts  in  that  valence  band  offsets  tend  to  be  smaller,  excitons  are  relatively 
strongly  bound  and  have  a  smaller  spatial  extension,  and  that  lasing  action  involves 
recombination  from  excitons  rather  than  band  states.  Systematic  investigation  of  these 
results  provides  some  guidance  as  to  the  magnitude  of  band  offsets  which  are  not  yet 
known  with  sufficient  accuracy.  The  small  valence  band  offsets  suggest  that  the  Te- 
based  superlattices  may  fall  into  a  weak  type  II  regime.  Variational  calculations  have 
been  performed  to  examine  the  behavior  of  excitons  in  this  regime.  These  calculations 
have  also  been  very  helpful  in  providing  confidence  that  the  simpler  binding  energy 
nonvariational  calculations  mentioned  above  are  quantitatively  accurate  for  type  I 
superlattices.  Because  of  the  considerable  lattice  mismatch  that  characterizes  these 
superlattices,  the  proper  incorporation  of  strain  effects  is  important. 

The  effects  of  applied  electric  and  magnetic  fields  on  superlattice  excitons  have 
also  been  systematically  investigated.  In  the  magnetic  field  case,  the  weak  field  limit, 
for  which  the  cyclotron  energy  is  less  than  the  mini-band  gap,  is  adequate  for  our 
purposes  since  the  weak  field  limit  extends  to  10-20  Tesla.  The  present  formalism 
suffices  for  the  interpretation  of  most  experimental  data  of  technical  interest.  The 
exciton  structure  is  observed  to  pass  smoothly  to  the  regime  in  which  a  description  in 
terms  of  Landau  levels  is  more  appropriate.  Theory  and  experiment  are  found  to  agree 
within  0.5  meV  for  the  energies  of  the  four  lowest  exciton  levels  in  fields  extending 
to  8  Tesla  in  two  different  InGaAs/GaAs  quantum  well  samples.  In  a  GaAs/GaAlAs 
superlattice  system  correspondingly  good  agreement  is  found  for  the  three  lowest  levels 
for  fields  up  to  20  Tesla.  The  exciton  binding  energy  increases  as  a  function  of  field, 
an  effect  which  is  most  pronounced  for  the  lowest  levels. 

Electric  fields  tend  to  localize  the  electron  and  hole  along  the  direction  of  the 
field  and  induce  the  formation  of  a  Stark  ladder  in  the  superlattice  at  fields  above 
about  104  V/cm.  The  present  formalism  is  applicable  to  fields  extending  up  to  10® 
V/cm.  For  higher  fields  the  independent  sub- band  approximation  breaks  down  due  to 
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interband  tunneling.  The  agreement  between  theory  and  experiment  is  excellent  for 
two  GaAs/GaAlAs  superlattices  that  have  been  examined  up  to  fields  of  120  kV/cm. 
The  binding  energy  of  the  strongest  exciton  is  enhanced  by  the  field  due  to  the  Stark 
localization  of  the  electron  and  hole  within  the  same  well. 

These  results  represent  part  of  P.M.  Young’s  Ph.D.  thesis  which  is  almost  com¬ 
pleted.  They  will  be  ready  for  publication  shortly. 

Observed  optical  structure  in  GaAs/GaAlAs  heterostructures  that  cannot  be 
explained  by  our  calculations  has  suggested  the  possibility  of  unintentionally  induced 
quantum  well  asymmetry  in  these  structures.  Such  asymmetries  would  break  the  usual 
optical  selection  rules  and  result  in  additional  structure.  Most  probably  these  asymme¬ 
tries  are  associated  with  unintentionally  introduced  electric  fields  created,  for  example, 
through  doping  imbalance  between  the  buffer  and  capping  layers  of  a  heterostructure. 
They  could  also  be  created  by  dangling  bonds  associated  with  the  intentional  miscut 
that  is  introduced  during  superlattice  growth.  Another  likely  possibility  is  that  com¬ 
positional  gradients  either  in  the  barriers,  the  wells  or  in  both  may  be  introduced  as  a 
result  of  the  MBE  growth  process.  For  example,  the  flux  of  the  A1  deviates  briefly  from 
its  nominal  value  when  shutters  are  open  or  closed  resulting  in  a  composition  gradient 
and  a  resulting  field.  All  observed  structure  can  be  explained  quantitatively  by  postu¬ 
lating  fields  of  the  order  of  104  V/cm.  The  results  have  been  accepted  for  publication 
in  Applied  Physics  Letters  [3].  The  simple  and  versatile  calculational  techniques  for 
calculating  optical  properties  that  have  been  developed  are  likely  to  be  useful  for  a 
variety  of  applications,  among  them,  the  characterization  of  actual  device  structures. 

References: 

1.  P.M.  Young,  P.M.  Hui  and  H.  Ehrenreich,  Pkys.  Rev.  B44,  12  969  (1991). 

2.  M.  Ziegler,  E.  Runge,  P.M.  Young  and  H.  Ehrenreich  (to  be  published). 

3.  P.M.  Young  and  H.  Ehrenreich,  Appl.  Pkys.  Lett,  (to  be  published). 
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B.  Response  end  Transit  Times  in  Quantum  Well  Structures 

The  electrical  properties  of  double  barrier  quantum  well  and  other  types  of  nanos¬ 
tructures  are  of  considerable  importance  because  of  their  potential  applicability  in 
computer  applications.  For  example  double  barrier  quantum  wells  (DBQW)  are  pro¬ 
totypical  circuit  components  in  cellular  automata  in  which  each  cell  is  required  to  have 
a  high  frequency  response  and  to  be  decoupled  from  all  but  its  nearest  neighbors.  Ef¬ 
fective  decoupling  can  be  achieved  if  wavefunction  coherence  between  adjacent  cells  is 
small.  The  question  addressed  in  the  present  work  concerns  the  relative  importance  of 
coherent  vs.  sequential  tunneling.  Phase  coherence  effects  are  not  present  in  the  latter 
process.  In  recently  published  work  [1,2]  the  response  of  a  bias  DBQW  structure  to  a 
small  ac  voltage  is  calculated  using  nonequilibrium  Green's  function  techniques.  The 
tunneling  is  shown  to  be  predominantly  sequential  in  InGaAs/InAlAs  structures  due 
to  only  well  alloy  scattering  which  is  described  in  the  coherent  potential  approxima¬ 
tion.  Of  particular  interest  is  the  fact  that  the  transit  or  dwell  time,  which  has  been 
measured  in  detail  by  Chemla  and  collaborators,  is  numerically  very  similar  to  the 
technologically  more  significant  ac  response  time,  which  is  far  more  difficult  to  mea¬ 
sure.  The  ac  response  time,  which  we  have  introduced,  provides  an  explicit  measure  of 
the  frequency  response  of  nanostructures  involving  tunneling  and  therefore  promises 
to  be  generally  useful.  Its  value  is  about  r  =  50  ps  for  the  structures  investigated. 

These  quantum  transport  calculations  have  been  extended  to  investigate  the  noise 
properties.  In  the  simplest  case  of  small  external  bias,  thermal  Johnson  noise  is  pre¬ 
dominant.  Shot  noise  is  predominant  for  larger  bias  when  the  quantum  well  structure 
as  a  whole  is  no  longer  in  thermal  equilibrium.  At  high  bias  voltages,  frequency  inde¬ 
pendent  shot  noise  results  for  systems  with  appreciably  asymmetric  barriers.  However, 
in  the  case  of  barriers  having  about  the  same  transparency,  the  noise  current  power 
spectrum  is  frequency  dependent  as  a  result  of  structure-induced  correlations  between 
subsequent  electron  tunneling  events.  This  frequency  dependence  is  characterized  by 
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the  same  ac  response  time  which  limits  the  high  speed  performance  of  such  devices. 
The  frequency  dependence  of  the  noise  current  power  spectrum  reflects  temporal  cor¬ 
relations  between  electrons  tunneling  into  and  out  of  the  well  [3]. 

Electron-electron  interactions  have  been  considered  in  connection  with  a  one¬ 
dimensional  model  describing  resonant  tunneling  through  a  quantum  dot  characterized 
by  a  ground  state  and  a  single  accessible  excited  state.  These  interactions  can  dominate 
the  device  behavior  of  nanostructures  operated  at  low  voltage.  In  order  to  understand 
the  associated  energy  and  time  scales  the  dynamic  response  is  evaluated  within  the 
context  of  a  Hubbard  Hamiltonian  using  nonequilibrium  Green’s  function  techniques. 

References: 

1.  E.  Runge  and  H.  Ehrenreich,  Phys.  Rev.  B45,  9145  (1992). 

2.  E.  Runge  and  H.  Ehrenreich,  Annals  of  Physics  (accepted  for  publication). 

3.  E.  Runge  (to  be  published). 

C.  Phonon  Assisted  Transport  Through  Tunneling  Structures 

The  goal  of  this  research  is  to  understand  the  origin  of  the  weak  replica  peak  ob¬ 
served  in  the  dc  I- V  curves  of  high  quality  DBQW  structures.  This  peak,  due  to  optical 
phonons,  appears  at  high  voltages  greater  in  magnitude  then  the  main  resonant  tun¬ 
neling  peak.  A  more  realistic  phonon  model  involving  confined  modes  with  dispersion 
is  introduced  to  describe  confined  and  interface  optical  modes.  The  present  calcula¬ 
tions,  based  on  the  same  kind  of  transport  calculations  described  in  Section  B  above, 
reproduce  the  observed  position  and  magnitude  of  the  replica  peak  and  confirm  its 
origin  as  arising  from  inelastic  scattering  of  electrons  due  to  quantum  well  vibrational, 
in  particular,  interfacial  modes.  The  use  of  the  virtual  crystal  bulk  phonon  spectrum 
in  these  calculations  would  spoil  the  good  agreement  with  experiment.  While  confined 
phonon  modes  have  been  previously  observed  Raman  scattering  experiments,  their 
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importance  in  transport  properties  has  been  previously  unappreciated  [1]. 

Reference: 

1.  C.  Grein,  E.  Runge  and  H.  Ehrenreich  (manuscript  in  preparation). 

D.  Auger  Processes  in  GalnSb/lnAs  and  HgCdTe  Infrared  Detectors  and  Quan¬ 
tum  Well  Lasers 

The  primary  focus  of  this  work  is  on  GalnSb/lnAs  superlattice  IR  detectors  which 
potentially  could  operate  at  a  longer  wavelength  with  reduction  in  tunneling  noise  and 
suppression  of  Auger  recombination.  Because  the  system  involves  III-V  rather  than 
II- VI  compounds  it  may  be  possible  to  fabricate  such  structures  more  reproducibly 
and  therefore  with  higher  yields.  The  goal  of  this  research  is  the  theoretical  evaluation 
of  “ideal”  device  characteristics.  This  involves  calculation  of  the  electronic  structure, 
optical  absorption  coefficients,  a  quantum  transport  theory  for  evaluating  electron  and 
hole  currents,  and  a  quantitative  assessment  of  electrical  noise. 

Substantial  funding  for  this  work  has  been  requested  from  DARPA/ONR.  Since 
many  of  its  results  are  also  applicable  to  quantum  well  lasers,  the  quantum  efficiency 
of  such  devices  will  also  be  explored  under  JSEP  auspices. 

As  may  be  gathered  from  the  research  efforts  described  above,  many  of  the 
techniques  for  undertaking  this  program  are  already  at  hand.  At  present  we  are  inves¬ 
tigating  nonradiative  Auger  processes  for  the  GalnSb/lnAs  superlattice  IR  detector 
on  the  basis  of  a  realistic  superlattice  band  structure.  We  have  obtained  preliminary 
results  and  identified  the  regions  of  K-space  predominantly  responsible  for  radiation¬ 
less  recombination.  The  computer  code  under  development  is  believed  to  be  the  most 
quantitative  tool  for  investigating  Auger  processes  presently  available. 

The  results  obtained  are  being  compared  to  HgCdTe  superlattice  and  alloy  IR 
detectors.  Preliminary  results  indicate  that  the  III-V  system,  which  is  being  investi- 
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gated  extensively  by  T.C.  McGill’s  group  at  Caltech  and  also  at  Hughes,  is  at  least 
comparable  to  HgCdTe  alloy  detectors.  The  preliminary  results  for  HgCdTe  superlat¬ 
tices  are  sufficiently  surprising  that  comments  concerning  them  will  be  deferred  to  a 
later  date  [1]. 

Reference: 

1.  C.  Grein,  P.M.  Young,  H.  Ehrenreich  and  T.C.  McGill  (submitted  to  MTC  Work¬ 
shop,  Boston,  1992). 
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1.2  Pressure  Dependence  of  Photoluminescence  and  Photoluminescence 

Excitation  in  Quantum  Wells  and  Superlattices.  J.H.  Burnett,  H.M. 

Cheong,  and  W.  Paxil,  Grant  N00014-89-J-1023;  Research  Unit  2. 

Using  hydrostatic  pressure  measurements,  we  have  obtained  information  on 
Group  Ill-Group  V  and  Group  II-Group  VI  quantum  well  and  superlattice  systems, 
which  have  significant  relevance  to  their  properties  at  atmospheric  pressure.  These  re¬ 
sults  are  firsts  of  their  kind,  and  we  believe  they  demonstrate  the  utility  of  hydrostatic 
pressure  measurements  for  understanding  the  properties  of  these  systems. 

A.  Group  II-Group  VI  Systems 

In  the  last  report,  we  reported  an  initial  measurement  of  the  photoluminescence 
(PL)  of  a  HgTe/CdTe  superlattice  (SL)  sample  at  pressures  up  to  15  kbar  at  80  K  in 
an  effort  to  test  the  current  theory  on  the  superlattice  band  structure  of  this  system. 
We  have  completed  this  project  with  an  extensive  set  of  measurements  on  another 
sample  at  pressures  up  to  30  kbar.  Three  major  PL  peaks  are  observed  in  the  100- 
700  meV  range.  Figure  1.1  shows  the  pressure  dependences  of  these  peaks.  The  most 
prominent  peak  at  ~  130  meV  has  a  pressure  coefficient  of  5  1  meV/kbar  as  reported 
in  the  previous  report.  The  calculated  energies  for  the  three  lowest  optical  transitions 
in  the  SL  at  zero  pressure  are  close  to  the  energies  of  these  peaks,  though  PL  from  the 
two  higher  energy  transitions  would  not  be  expected  to  be  observable  at  80  K.  The 
higher  energy  PL  peaks,  which  have  similar  pressure  coefficients  (1-2  meV/kbar)  as 
the  lowest  energy  one,  are  not  from  the  barrier  material  (CdTe)  or  from  the  substrate 
(CdZnTe),  for  these  peaks  were  absent  in  the  PL  spectra  of  these  materials. 

Mr.  Paul  Young  of  Prof.  Ehrenreich’s  group  has  carried  out  calculations  of  the 
pressure  dependence  of  the  SL  bands  within  tin  envelope  function  approximation 
(EFA)  based  on  the  eight  band  Kane  model,  using  a  valence  band-offset  (A)  value 
of  350  meV.  Calculations  indicate  that  the  SL  band  gap  increases  with  pressure  at  a 
rate  of  ~  6.5  meV/kbar  near  zero  pressure  and  this  rate  increases  to  ~  8  meV/kbar 
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Figure  1.1.  Pressure  dependences  of  the  three  moat  prominent  photbluminescence  peaks  from 
a  HgTe/CdTe  sample.  The  inset  shows  a  typical  luminescence  spectrum,  taken  at  20.4  kbar. 
All  three  peaks  have  pressure  coefficients  of  1-2  meV/kbar. 
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when  the  pressure  reaches  20  kbar.  Varying  the  input  parameters  within  reasonable 
ranges  changes  the  pressure  dependences  by  less  than  10%.  A  small  A  (40  meV)  de¬ 
creases  the  pressure  dependence  by  less  than  10%.  The  next  lowest  energy  transitions 
have  similar  pressure  dependences. 

The  pressure  dependence  of  the  SL  band  gap  of  >  6  meV/kbar,  calculated  in  the 
EFA,  is  far  outside  the  error  bars  of  the  measured  pressure  dependence  of  the  main 
peak,  <  1  meV/kbar.  The  disagreement  between  theory  and  experiment  suggests  at 
least  one  of  three  possibilities.  First,  the  EFA  may  be  inappropriate  to  determine 
the  band  gap  and  its  pressure  dependence.  This  approximation  has  been  effective 
in  determining  the  band  gap  and  its  pressure  dependence  for  GaAs/AlrGai_xAs  het¬ 
erostructures,  but  it  may  not  be  valid  in  describing  HgTe/CdTe  SL’s  due  to  the  unusual 
band  structure  of  HgTe. 

Second,  the  assignment  of  the  strong  PL  peak  at  atmospheric  pressure  to  the  fun¬ 
damental  SL  band  gap  may  be  incorrect.  Because  some  of  the  SL  bands  are  so-called 
“quasi-interface  bands,”  which  have  large  probability  densities  at  the  HgTe/CdTe  in¬ 
terface  where  a  laxge  number  of  interface  defects  also  exist,  the  fundamental  PL  may 
not  be  strong  enough  to  be  seen  and  the  observed  PL  may  be  defect  luminescence. 
This  is  supported  by  the  fact  that  the  PL  peak  and  the  absorption  edge  have  differ¬ 
ent  temperature  dependences.  If  this  is  the  case,  the  PL  of  HgTe/CdTe  may  have 
been  generally  misidentified  in  the  previous  studies.  Absorption  measurements  as  a 
function  of  pressure  would  be  a  decisive  test  of  this  possibility.  However,  it  should 
be  noted  that  generally,  defect  luminescence  peaks  axe  lower  in  energy  than  the  band 
gap  luminescence,  and  pressure  dependences  of  defect  luminescence  peak  energies  are 
usually  similar  to  that  of  the  band  gap. 

Third,  one  of  the  assumptions  ancillary  to  the  use  of  the  EFA  may  be  incorrect. 
This  could  be  the  assumption  that  it  is  proper  to  ignore  the  possibility  that  a  defected 
atomic  arrangement  at  the  HgTd/ CdTe  interface  is  determining  one  or  other  of  the 
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band  edge  states.  It  could  also  be  the  assumption  that  A  at  ambient  pressure  is 
350  meV,  although  this  assumption  is  currently  thought  to  be  firmly  based.  And 
finally,  the  assumption  of  pressure  independence  of  A  could  be  wrong.  Even  though 
this  assumption  is  generally  believed  to  be  valid,  with  some  theoretical  support,  there 
has  been  no  experimental  proof,  and  there  has  been  a  suggestion  of  a  temperature 
dependence  of  A.  Our  calculations  of  the  effect  due  to  the  pressure  dependence  of 

A,  however,  indicate  that  to  reduce  the  pressure  dependence  of  the  gap  to  4  and 
1  meV/kbar,  A  has  to  increase  by  15  and  30  meV/kbar,  respectively,  which  seems 
unreasonably  high. 

We  plan  to  measure  the  absorption  spectra  as  a  function  of  pressure.  This  would 
eliminate  the  question  of  whether  the  PL  is  from  the  transitions  between  the  SL  bands. 
If  the  absorption  has  the  same  pressure  dependence  as  the  PL,  then  the  validity  of  the 
EFA  and/or  its  ancillary  assumptions  will  have  to  be  reconsidered. 

B.  Group  Ill-Group  V  Systems 

One  of  the  most  important  parameters  for  understanding  and  designing  het¬ 
erostructure  systems  is  the  relative  energy  of  the  conduction  band  (CB)  edges 
at  the  junctions  between  the  materials,  the  conduction  band  offset.  For  the 
GaAs/AlIGa1_xAs  system,  probably  the  most  reliable  value  determined  for  this  pa¬ 
rameter,  QcB&Eg  =  0.69 &Eg  where  AEg  is  the  band  gap  difference,  was  determined 
by  hydrostatic  pressure  measurements,  by  measuring  the  pressure  at  which  the  X  band 
of  the  AljGa^-xAs  layer  crosses  the  CB  minimum  of  the  GaAs  layer,  resulting  in  a 
separation  of  the  electron  and  holes,  and  thus  a  dramatic  reduction  of  the  photolumi¬ 
nescence  intensity  [1].  However,  this  value  depends  on  one  crucial  assumption,  namely 
that  the  band  offset  is  itself  independent  of  pressure.  Until  now  there  has  been  no  di¬ 
rect  evidence  for  this  assumption.  We  have  now  completed  preliminary  measurements 
on  an  experiment  to  put  this  assumption  to  direct  test. 
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To  determine  the  pressure  dependence  of  the  conduction  band  offset,  we  have 
measured  the  harmonic  oscillator  energy  levels  of  a  1000  A-wide  Al*Gai_r  As  parabolic 
quantum  well,  using  photoluminescence  excitation  (PLE)  spectroscopy,  as  a  function 
of  hydrostatic  pressure  up  to  35  kbar.  The  sample  was  grown  and  characterized  by 
Professor  A.  Gossard  at  the  University  of  California,  at  Santa  Barbara,  CA.  A  typical 
PLE  spectrum  of  the  parabolic  well,  here  taken  at  7.1  kbar,  is  shown  in  the  inset  of 
Figure  1.2,  with  the  principal  transition  labeled.  The  harmonic  oscillator  spacing  is 
obtained  from  differences  in  the  transition  energies.  In  a  simple  model  the  harmonic 
oscillator  spacing  is  given  by: 


hu  = 


'8  QcB&Eg 
L\m* 


where  Lx  is  the  width  of  the  parabolic  well,  &.Eg  is  the  difference  between  the  band 
gap  energies  at  the  edge  and  the  center  of  the  well,  m*  is  the  average  electron  effective 
mass,  and  Qcb  is  the  fraction  of  the  band  gap  difference  in  the  conduction  band.  Since 
hu  depends  directly  on  Qcb,  we  can  directly  determine  the  pressure  dependence  of 
Qcb  by  measuring  hu  as  a,  function  of  pressure.  Our  measurement  of  the  pressure 
dependence  of  hu  is  shown  in  Figure  1.2.  Also  included  are  calculations  for  the  pressure 
dependence  of  hu ,  including  the  pressure  effects  on  I,  and  m*,  for  several  assumed 
pressure  dependences  of  the  band  offset  Qcb,  dQcB/dP  (1/kbar).  (In  the  top  curve 
the  pressure  dependence  of  the  effective  mass  is  neglected.)  Our  measurements  indicate 
that  Qcb  changes  less  than  7%  for  pressures  from  0  to  35  kbar.  If  these  results  are 
confirmed  on  further  measurements,  it  supports  the  validity  of  the  value  for  the  band 
offset  determined  by  pressure  measurements  [1].  It  also  has  significance  on  the  role  of 
strain  on  the  band  offsets  in  strained  layer  superlattices. 
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Figure  1.2.  Preaaure  dependence  of  the  conduction  band  harmonic  oedllator  spacing  of  a 
1000  A  Al*Gai_x  As  parabolic  well,  determined  by  measuring  the  pressure  dependence  of  the 
energy  difference  between  the  photoiuminescence  peaks  H1-E3  and  Hl-El  shown  in  the  inset. 
The  lines  show  the  calculated  pressure  dependences,  including  the  pressure  effects  on  L,  and 
mj,  for  several  assumed  pressure  dependences  of  the  band  offset  QcBi  4Qcb/4P  (1/kbar). 
(In  the  top  curve  the  pressure  dependence  of  the  effective  mass  is  neglected). 
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I.S  X-ray  Surface  Characterization.  I.M.  Tidswell,  T.  Rabedeau,  and  P.S.  Per* 
shan,  Grants  N00014-89-J-1023,  NSF  DMR-89-20490,  NSF  DMR  91-13782; 

Research  Unit  3. 

The  results  on  the  X-ray  scattering  study  of  the  SiO2/Si(001)  interface  that  was 
formed  by  dry  oxidation  of  the  Si(001)  were  described  in  the  report  for  the  previous 
period  and  the  publications  resulting  from  that  work  have  now  appeared  in  print 
[1,  2].  That  project  had  been  stimulated  by  published  disagreement  over  whether 
or  not  the  SiC>2/Si(001)  interface  consisted  of  an  epitaxial  layer  of  ordered  oxide  [3- 
8].  Our  results  resolved  some  of  this  controversy  by  demonstrating  that  for  a  room 
temperature  oxidation  process  the  epitaxially  ordered  fraction  of  Si02  covers  no  more 
than  —  8%  of  the  interface,  the  remaining  being  amorphous,  and  even  this  small 
fraction  of  epitaxially  ordered  Si02  is  unstable  against  either  exposure  to  ambient  air 
(i.e.,  wet  O2)  or  1  hour  post-oxidation  anneal  at  550°C.  We  conjectured  that  the  low 
initial  fractional  coverage  and  the  instability  against  humidity  or  temperature  might 
have  been  related  to  the  inherent  strain  induced  by  the  lattice  mismatch  between 
the  Si(001)  and  the  SiC>2  [9-12],  and  we  proposed  studying  the  oxidation  of  Si  layers 
that  are  epitaxially  deposited  on  substrates  made  of  alloys  of  SiIGej_,/Ge(001).  By 
varying  the  concentration  over  the  range  1  >  1  >  0  tensile  strains  in  the  Si  lattice 
can  be  varied  from  0  to  >  4%.  Since  application  of  an  external  compressive  stress  is 
known  to  reduce  the  internal  stress  on  the  dimer,  we  expect  that  the  dimer  interfacial 
structure  will  be  less  stable  when  deposited  on  the  Si*Geio*/Ge  substrate  than  on  the 
Si  [13]. 

Although  sample  preparation  for  this  project  was  delayed  some  months  as  a 
result  of  necessary  maintenance  for  the  particular  AT&T  MBE  machine  that  this 
program  was  designed  around,  the  machine  was  recommissioned  at  the  end  of  1991 
and  a  good  Ge(001)  surface  was  prepared  and  evaluated  using  the  Harvard  Rotating 
Anode.  To  form  that  sample  a  single  crystal  (miscut  <  0.5°)  was  epitaxially  coated 
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Figure  1.3.  The  ratio  of  the  measured  X-ray  reflectivity,  R(Qa)  for  the  Ge[001]/Si/SiOj 
described  in  the  text,  to  the  theoretical  Fresnel  reflectivity,  Rf(QM)  for  a  flat  surface  with  the 
same  average  electron  density  as  crystalline  Ge.  The  minima  at  QM  2U  0.28 A-1  and  0.88 A-1 
correspond  to  destructive  interference  between  reflections  from  the  2!  llA  thick  Si/SiOj  layer 
(e.g.,  Qad  =  w  and  3*j.  Note  that  the  interference  minima  at  QM  9!  1.45A-1,  Qai  —  Sw, 
corresponding  to  the  anti-Bragg  (002)  point  of  Ge  is  missing  as  expected. 


with  approximately  2000A  of  Ge  at  400°C,  following  which  it  was  annealed  at  500°C 
for  10  minutes;  four  monolayers  of  Si  were  then  epitaxially  deposited  between  250- 
350°C.  The  sample  was  exposed  to  dry  O2  at  room  temperature  for  approximately  20 
hours. 

The  results  of  specular  reflectivity  measurements  that  were  made  in  March  1992 
on  beam  line  X20  at  the  National  Synchrotron  Light  Source  (NSLS)  are  shown  in 
Figure  1.3. 

A  second  beam  allotment  on  X20a  at  NSLS  was  assigned  to  this  project  for  the 
period  from  May  14th  through  May  31st.  Since  this  report  is  being  prepared  on  June 
1  we  present  only  the  data  with  a  minimum  of  analysis.  We  have  used  the  splitting 
between  four  inequivalent  (1,1, ±1)  and  (1,3,±1)  truncation  rods  near  to  the  [001]  plane 
to  measure  the  angle  of  the  miscut  between  the  surface  normal  and  the  [001]  axis  of 
the  same  Ge/Si/Si02  sample  that  was  studied  in  March;  it  is  approximately  0.05°  and 
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near  to  the  [100]  direction.  In  contrast  to  the  SifOOlJ/SiOa  surface  that  we  previously 
studied  we  do  not  observe  long-range  order  in  the  surface  steps  [2].  Presumably  this 
is  because  in  the  earlier  experiment  the  miscut  for  the  Si  crystal  was  close  to  the  [110] 
direction  that  is  favored  by  surface  steps  on  the  [001]  face  of  Si  and  Ge  [14-19].  With 
the  miscut  pointing  towards  the  [100]  direction  the  surface  steps  are  surely  faceted  with 
individual  segments  along  the  [1,±1,0]  directions  and  long-range  order  is  less  likely. 

We  have  measured  eight  different  fractional  order  surface  peaks  and  demon¬ 
strated  that  the  anisotropy  of  the  surface  correlation  lengths  are  significantly  shorter 
along  the  miscut  direction  than  in  the  perpendicular  direction.  Figure  1.4  illustrates 
the  raw  data  for  the  HWHM  of  transverse  (e.g.,  rocking)  scan  through  these  eight 
fractional  order  peaks  as  a  function  of  the  angle  between  the  direction  of  the  scan  and 
the  [100]  axis.  The  broken  line  was  hastily  drawn  to  illustrate  the  systematics  of  the 
data.  It  represents  0.003  x  |  cos[^  -  <j>]\  where  'P  is  the  direction  of  the  transverse  scan 
through  a  surface  peak,  e.g.,  for  a  transverse  scan  parallel  to  the  surface  and  through 
the  (1/2,  3/2,  .03)  point,  the  angle  $  =  tan-1(l/-3)  =  108.43°.  The  angle  of  the 
miscut  away  from  the  [100]  direction  was  calculated  from  the  observed  splitting  to  be 
</•  =  5°. 

The  Qz  dependence  of  the  (1,1, Qz)  truncation  rod  and  the  (1/2, 1/2, Qx)  frac¬ 
tional  order  rod  have  been  measured.  Interference  oscillations  on  the  (1,1, Q,)  rods 
and  the  absence  of  any  structure  on  the  (1/2, 1/2 ,QX)  suggest  that  the  fractional  order 
peaks,  which  we  presume  to  be  due  to  the  1x2  dimer  previously  seen  on  the  surface 
oxidized  surface  of  Si(001)/SiC>2  [20],  are  only  one  monolayer  thick  while  the  com¬ 
mensurate  (lxl)  order  of  the  epitaxial  SiC>2  may  be  four  layers  thick.  Quantitative 
analysis  will  be  confirm  some  details  of  the  initial  interpretation  of  the  raw  data;  there 
are  a  number  of  further  experiments  that  will  be  necessary  in  order  to  obtain  a  firm 
picture  of  this  surface. 
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Figure  1.4.  Illustration  of  the  dependence  of  the  observed  HWHM  for  transverse  GID  scans 
through  eight  different  fractional  order  peaks.  The  direction  that  a  transverse  scan  through  a 
peak  (*,*, 0.033)  makes  with  respect  to  the  (100)  direction  is  given  by  ♦  =  ta n~l{h/—k).  The 
broken  line,  which  is  simply  a  guide  for  the  eye  is  given  by  0.003  x  |coe[¥  —  ^]|  where  4  —  5°  is 
the  angle  between  the  projection  of  the  miscut  direction  on  the  surface  and  the  [001]  direction. 

The  most  compelling  X-ray  measurements  that  need  to  be  made  on  the  present 
sample  include  characterization  of  other  fractional  order  rods.  It  is  not  proven  that 
the  surface  domains  are  (1  x  2)  rather  than  (2  x  2)  [20];  however,  on  the  basis  of  that 
assumption  we  would  like  to  see  if  the  two  inequivalent  (1x2)  domains  are  equally 
occupied.  To  that  end,  and  also  to  get  further  information  about  the  structure  factor 
of  the  (1x2)  surface  lattice,  we  will  also  measure  scattering  along  the  Qz  direction  for 
other  fractional  order  rods.  We  would  also  like  to  compare  the  specular  reflection  with 
the  miscut  oriented  perpendicular  to  the  plane  of  incidence  with  the  data  in  Figure  1.3 
which  was  taken  with  the  miscut  in  the  plane  of  incidence.  This  will  give  information 
on  the  degree  of  microroughness  within  the  individual  terraces  in  the  stepped  surface. 
The  most  important  measurement  will  be  to  tune  the  beam  line  to  an  energy  near 
to  the  peak  in  the  anomalous  dispersion  of  the  Ge  scattering  factor,  A  2  1.13A  in 
order  to  compare  with  scattering  observed  at  A  =  1.41  A  that  is  currently  being  used. 
Since  the  near  resonant  energy  will  enhance  the  scattering  amplitude  of  the  Ge  by 
nearly  a  factor  of  10  relative  to  that  of  the  Si  and  0  this  experiment  would  enable 
us  to  produce  independent  models  for  the  Si  and  Ge  surface  distributions.  We  will 
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also  explore  possibilities  of  using  other  techniques,  such  as  TEM,  STM,  SEXAFS, 
and  photoelectron  spectroscopy  to  obtain  information  that  would  complement  our 
structural  results. 


At  some  point  we  will  prepare  Si  layers  on  alloy  substrates  such  as  SirGei_*  and 
make  comparable  studies  in  order  to  systematically  vary  the  effect  of  strain  on  the 
SiOj  surface  layers. 
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1.4  Superconducting  Josephson  Junction  Arrays.  M.S.  Rz  chow  ski,  L.L. 
Sohn,  T.S.  Tighe,  and  M.  Tinkham,  Grants  N00014-89-J- 1023,  N00014-89- 

J-1565,  and  DMR-89-12927;  Research  Unit  4. 

After  winding  up  the  work  on  giant  Shapiro  steps  stemming  from  coherent  fluxon 
motion  induced  by  rf-drive  currents  in  conventional  SNS  arrays  of  various  geometries, 
we  have  shifted  our  focus  to  a  novel  array  geometry  which  results  in  long-range  interac¬ 
tions.  Rather  than  consisting  of  superconducting  islands  which  are  Josephson- coupled 
to  nearest  neighbors,  these  new  arrays  consist  of  two  orthogonal  sets  of  N  parallel 
superconducting  wires  which  are  coupled  by  a  Josephson  junction  at  every  point  of 
crossing.  (See  Figure  1.5.)  As  a  result  of  this  unique  geometry,  all  the  wires  in  the  ar¬ 
ray  are  Josephson- coupled  to  each  other  as  nearest-  or  next-nearest  neighbors.  Instead 
of  being  coupled  to  4  or  6  neighbors,  they  axe  coupled  to  N  (~  1000)  neighbors,  and  a 
mean-field  description  might  be  expected  to  be  applicable.  Indeed,  this  approach  was 
used  by  Vinokur  et  al.  [1],  who  first  drew  attention  to  the  interesting  properties  of 
such  arrays,  particularly  in  the  presence  of  a  magnetic  field  and  if  the  spacing  of  the 
wires  was  disordered. 

Our  approach  to  this  system  has  combined  theory  and  experiment.  On  the 
theory  side,  we  worked  out  the  phase  transition  of  the  system  to  a  state  showing  long- 
range  phase  coherence  [2].  We  did  this  both  by  a  mean-field  approximation  and  by 
Monte  Carlo  simulations,  for  regular  and  disordered  arrays,  and  both  in  and  out  of  a 
magnetic  field.  Our  work  confirmed  that,  in  zero-field,  the  phase  transition  to  a  phase- 
coherent  state  should  occur  at  a  temperature  Tc  which  scales  with  N,  the  number  of 
wires  coupled  as  nearest-neighbors  to  any  given  wire.  In  the  presence  of  a  field,  Tc  is 
depressed  by  an  amount  depending  on  the  field  strength,  and  upon  whether  the  arrays 
are  ordered  or  disordered. 

On  the  experimental  side  [3],  we  developed  a  technique  to  fabricate  such  arrays, 
typically  with  N  —  1000  parallel  Nb  thin-film  wires  in  each  set,  with  a  2.5  pm  line  width 


24 


Superconducting  Wires 


Josephson 

Junction 


(a)  (b) 

Figure  1.5.  Sections  of  an  (a)  ordered  and  (b)  disordered  array.  The  black  lines  are  the 

superconducting  wires.  Each  wire  is  coupled  to  every  orthogonal  wire  by  a  Josephson  junction. 

and  5  nm  spacing,  and  with  AlOx  barrier  tunnel  junctions  at  each  point  of  crossing. 
These  individual  junctions  typically  had  normal  state  resistances  rn  of  ~  670  Q  and 
critical  currents  ic  of  ~  0.5  pA.  When  we  measured  the  collective  properties  of  the 
arrays,  however,  we  found  that  they  were  far  from  those  predicted  by  the  idealized 
model.  For  example,  the  critical  current  between  one  wire  and  another  was  not  Nic 
as  would  be  expected  naively  for  N  (=1000)  junctions  in  parallel,  but  Ntffic,  where 
Neff  ~  20  <<  1000.  Similarly,  the  Tc  at  which  long-range  phase  coherence  sets  in  also 
acts  as  if  an  even  smaller  number  Ntff  of  wires  were  relevant.  We  have  traced  these 
major  discrepancies  from  standard  theory  to  the  fact  that  a  phase-gradient  is  set  up 
along  the  drive  wire  by  the  current  itself,  in  addition  to  that  caused  by  any  externally 
applied  field.  This  sets  up  a  “soliton”  structure  in  the  current  pattern,  which  limits  jVeff 
to  a  number  that  scales  with  the  square  root  of  the  ratio  of  the  inductance  of  the  wire  in 
one  unit  cell  to  the  Josephson  inductance  of  a  junction.  We  have  developed  a  detailed 
model  which  accounts  well  for  the  reduction  in  the  critical  current  of  the  array  below 
the  naive  estimate,  and  accounts  at  least  qualitatively  for  the  other  discrepancies.  This 
work  points  the  way  to  a  second  generation  set  of  experiments,  with  improved  samples 


designed  taking  account  of  these  inductive  effects.  It  also  points  to  the  need  for  a  more 
comprehensive  theory  of  this  novel  and  interesting  system  in  the  dynamic  state  above 
Ic,  where  surprising  hysteretic  and  nonlinear  effects  are  observed  experimentally. 
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Rev.  B,  submitted. 


1.5  Single-Electron  Tunneling  Devices  A.  Hanna,  M.  Tuominen,  J.  Hergen- 
rother,  and  M.  Tinkham,  Grants  N00014-89-J-1023,  N00014-89-J-1565,  and 
DMR-89-12927;  Research  Unit  4. 

These  devices  are  based  on  the  use  of  mesoscopic  structures  which  are  sufficiently 
small  that  their  capacitance  C  is  less  than  10-15  F,  so  that  the  charging  energy  e2/2 C 
of  a  single  electron  is  the  dominant  energy  in  the  problem.  We  are  pursuing  two 
approaches  to  this  regime:  one  based  on  the  STM,  and  the  other  based  on  thin  film 
nanofabrication. 

The  first  approach  to  succeed  was  the  STM-based  geometry.  In  this  there  is 
a  gold  grain  of  ~  50A  diameter,  into  which  current  tunnels  from  the  STM  tip,  and 
then  tunnels  out  to  the  Nb/NbOx  substrate.  With  this  device,  capacitances  as  small 
as  10-18  F  zue  readily  obtained.  This  is  small  enough  that  the  device  shows  little 
thermal  rounding  of  its  I-V  characteristics  at  4K,  and  shows  observable  Coulomb 
blockade  effects  even  at  room  temperature.  Our  first  successes  [1]  with  this  technique 
of  forming  a  “single-electron  transistor”  were  reported  in  the  previous  Annual  Report, 
confirming  that  a  charge  shift  of  only  e/2  is  sufficient  to  take  the  device  from  zero  to 
maximum  conductance.  Our  major  advance  during  the  past  year  was  to  demonstrate 


26 


for  the  first  time  [2]  the  observability  of  elastic  q-MQT,  i.e.,  the  so-called  macroscopic 
quantum  tunneling  of  charge  by  a  coherent  double  tunneling  process  which  leaves  no 
excitation  energy  behind  in  the  grain  between  the  two  tunnel  links.  This  process  yields 
a  tiny  linear  conductance  in  the  Coulomb  blockade  region,  ~  1000  times  smaller  than 
outside  of  the  Coulomb  blockade,  but  observable  at  the  pA  level  of  sensitivity.  Since 
the  strength  of  this  process  depends  inversely  on  the  volume  of  the  grain,  the  tiny 
grain  used  in  this  configuration  is  essential  to  the  observability  of  this  effect,  which 
has  never  been  seen  in  a  fabricated  thin  film  device. 

In  parallel  with  the  STM  work,  we  have  been  developing  techniques  for  nanofab¬ 
rication  of  thin-film  realizations  of  the  single-electron  tunneling  transistor  concept. 
We  use  an  AI-AI2O3-AI  structure,  with  600A  line  width,  yielding  junctions  with  ca¬ 
pacitance  2  x  10-16  F.  This  is  low  enough  to  allow  operation  at  temperatures  which 
can  be  reached  without  use  of  a  dilution  refrigerator,  although  we  do  use  one  to  ex¬ 
plore  the  limiting  behavior  at  very  low  temperatures.  Our  major  focus  in  this  work 
is  to  emphasize  the  unique  features  which  arise  when  the  A1  is  in  the  superconduct¬ 
ing  state,  ;n  which  there  is  a  competition  between  the  superconducting  energy  gap 
and  the  Coulomb  charging  energy  gap.  This  causes  the  current  through  the  device  to 
show  a  complex  periodic  variation  with  the  bias  potential  on  the  gate  electrode.  In  all 
previous  work  this  variation  has  shown  an  e-periodicity,  repeating  each  time  the  gate 
voltage  is  changed  sufficiently  to  pull  smother  electron  onto  the  island.  However,  in  one 
particularly  good  sample,  we  have  found  [3]  a  2e-periodicity  for  all  T  <  300  mK.  The 
crossover  in  period  is  shown  in  Figure  1.6.  This  2e  period  reflects  the  dominance  of 
Cooper  pairs  in  the  superconducting  state,  but  in  less  perfect  samples  single-electron 
tunneling  events  obscure  its  presence.  We  can  account  theoretically  for  our  data  by 
noting  that,  at  low  temperatures,  the  superconducting  pairing  energy  favors  states 
with  an  even  number  of  electrons  on  the  island.  The  cut-off  at  300  mK  reflects  the 
fact  that  at  this  temperature,  statistically  speaking,  a  single  quasiparticle  is  thermally 
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Figure  1.6.  Dependence  of  device  current  on  gate  voltage  at  biaa  voltage  of  100  pV  and 
varioua  temperatures.  Other  data  not  shown  here  show  that  the  2 e- periodicity  vanishes  near 
300  mK. 

excited  above  the  gap  in  the  entire  island  electrode.  We  believe  this  is  the  first  time 
that  an  experimental  distinction  has  been  made  between  an  odd  and  an  even  number 
of  electrons  on  a  fabricated  meted  structure. 
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d.  Books  (and  sections  thereof)  Published 

M.  Tinkham,  “Josephson  Effect  in  Low-Capacitance  Tunnel  Junctions,”  Series  of 
Invited  Lectures  at  NATO  ASI  on  Single  Charge  Tunneling,  Les  Houches,  March 
1991;  chapter  in  Single  Charge  Tunneling,  Grabert  and  Devoret,  eds.,  Plenum, 
NY,  1992,  pp.  139-166. 

M.  Tinkham,  “Conference  Summary,”  Los  Alamos  Symposium  1991,  published  in 
Phenomenology  and  Application  of  High-Temperature  Superconductors,  Addison- 
Wesley,  1992,  pp.  499-514. 

M.  Tinkham,  “Josephson  Effects  in  Weak  Links,”  article  in  Concise  Encyclopedia 
of  Magnetic  and  Superconducting  Materials,  edited  by  J.  Evetts,  Pergamon  Press, 
1992,  pp.  203-207. 

h.  Contributed  Presentations  at  Topical  or  Scientific/Technical  Society 
Conferences 

M.S.  Rzchowski,  L.L.  Sohn,  and  M.  Tinkham,  “Josephson  Junction  Arrays  with 
Long-Range  Interaction,”  March  1992  Meeting  of  the  American  Physical  Society. 

L.L.  Sohn,  M.S.  Rzchowski,  M.T.  Tuominen,  J.U.  Free,  and  M.  Tinkham,  “Prop¬ 
erties  of  Josephson  Junction  Arrays  with  Long-Range  Interaction,”  March  1992 
Meeting  of  the  American  Physical  Society. 

T.S.  Tighe,  M.T.  Tuominen,  J.M.  Hergenrother,  and  M.  Tinkham,  “Charging 
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A.E.  Hanna  and  M.  Tinkham,  “Single  Electron  Dynamics  in  a  Two  Junction 
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Dr.  M.S.  Rzchowski  and  Mr.  N.  Anand. 

Note:  Five  other  students  share  JSEP  facilities,  but  have  no  JSEP  salary.  Also, 
all  listed  publications  acknowledge  shared  support  from  the  other  grants  enumer¬ 
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1.6  Neural  Networks.  F.R.  Waugh,  R.  Seshadri,  C.  Marcus,  and  R.M.  Wester- 

velt,  Grant  N00014-89-J-1023;  Research  Unit  5. 

For  a  number  of  years  we  have  focussed  our  efforts  in  the  Joint  Services  Electron¬ 
ics  Program  on  electronic  neural  networks,  and  on  the  dynamics  of  magnetic  domain 
patterns  in  garnet  films.  This  research  has  been  quite  successful  as  described  below  and 
in  previous  annual  reports.  Beginning  this  year,  we  intend  to  shift  our  JSEP-supported 
efforts  toward  the  physics  of  electronic  nanostructures.  With  recent  improvements  in 
electron  beam  lithography,  we  now  have  state-of-the-art  facilities  in  the  McKay  labora¬ 
tory  for  the  production  and  characterization  of  nanostructures.  Professor  Westervelt’s 
group  uses  these  facilities  to  make  and  study  GaAs/AlGaAs  nanostructures  in  collab¬ 
oration  with  Prof.  A.C.  Gossard  at  the  University  of  California,  Santa  Barbara,  while 
Professor  Tinkham’s  group  is  expert  in  the  fabrication  and  study  of  metallic  and  su¬ 
perconducting  nanostructures.  The  presence  of  both  groups  in  the  McKay  laboratory 
offers  excellent  opportunities  for  new  collaborative  work,  because  many  of  the  physical 
issues  are  similar.  For  example,  as  described  below  we  are  currently  developing  hy¬ 
brid  semiconductor /superconductor  samples  in  which  a  Josephson  junction  is  used  to 
illuminate  a  semiconducting  quantum  point  contact  with  submillimeter  wavelength  ra¬ 
diation.  If  successful,  this  technique  will  permit  spectroscopy  of  single  nanostructures 
in  this  difficult  frequency  range. 

A.  Competitive  Neural  Networks  for  Vision 

In  the  past  we  have  developed  general  techniques  to  stabilize  fully  interconnected 
electronic  neural  networks  by  clocking,  and  have  conducted  a  systematic  study  of  their 
dynamics  and  stability.  In  the  past  year  we  have  applied  these  techniques  by  developing 
a  new  type  of  neural  net  architecture  for  problems  in  vision.  The  general  approach  is 
similar  to  the  neocognitron  developed  by  Fukushima.  The  network  is  subdivided  into 
identical  cells  which  are  programmed  to  recognize  elementary  image  features.  These 
features  are  then  combined  via  inter-cell  connections  into  more  complex  objects.  In 
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competitive  networks,  we  represent  each  elementary  feature  by  a  single  neuron  in  a 
cell  and  implement  a  variation  of  the  winner-take- all  circuit  by  requiring  that  the 
summed  output  of  all  neurons  within  a  given  cell  is  a  constant.  This  architecture 
requires  many  fewer  interconnections  than  our  previous  approach  in  which  each  cell 
was  a  fully  interconnected  network  of  pixels,  and  is  well  suited  to  implementation 
in  hardware.  Assembly  of  the  elementary  features  favored  by  individual  cells  into  a 
single  recognized  object  is  done  via  inter-cell  connections  which  can  either  be  all-to-all 
or  limited  to  nearest  neighbors.  We  have  conducted  a  systematic  analytical  study 
of  the  stability  of  clocked  neural  nets  of  this  type,  and  demonstrated  the  concept  in 
numerical  experiments  [1]. 

B.  Statistical  Mechanics  of  Magnetic  Bubble  Arrays 

In  the  past  year  we  have  completed  a  study  of  the  dynamics  and  statistical 
mechanics  of  magnetic  bubble  arrays  in  thin  garnet  films  [2-5].  A  magnetic  bubble 
array  is  an  experimentally-accessible  example  of  a  two-dimensional  fluid  similar  in 
some  respects  to  vortex  arrays  in  type-II  superconductors.  We  use  observations  of  the 
dynamics  of  bubble  arrays  to  test  theories  of  vortex  dynamics  in  collaboration  with 
Professors  David  Nelson  and  Daniel  Fisher  at  Harvard.  We  have  shown  that  bubble 
arrays  undergo  a  hexatic-to-liquid  melting  transition  [2],  and  have  systematically  stud¬ 
ied  their  statistical  mechanics  [3-4].  Most  recently  we  have  studied  shear  flow  past  a 
planar  interface,  in  order  to  test  the  predictions  of  theories  of  flux  pinning  in  high  Tc 
superconductors  [5].  Flux  flow  is  important  in  high  Tc  materials,  because  it  produces 
power  dissipation  and  destroys  the  advantage  of  superconductivity.  This  work  was 
supported  primarily  by  a  separate  grant  from  the  Office  of  Naval  Research,  and  we 
will  shift  the  support  of  future  work  in  this  area  entirely  to  the  ONR. 

C.  Semiconductor  Nanostructure* 

In  the  past  year  we  have  developed  and  used  new  facilities  for  the  production 
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and  characterization  of  semiconductor  nanostructures,  with  partial  JSEP  support. 
Electron  beam  lithography  has  been  used  to  make  submicron  quantum  boxes  with 
quantum  point  contacts  which  we  have  used  to  study  chaotic  scattering  of  electrons  [6]. 
We  find  that  chaotic  scattering  of  electrons  from  the  walls  of  a  quantum  box  produces 
features  in  the  low  temperature  magnetoresistance  analogous  to  weak  localization  and 
universal  conductance  fluctuations.  The  difference  is  that  scattering  in  our  structures 
is  primarily  from  the  walls  of  the  box  and  depends  on  its  shape,  while  localization  and 
universal  conductance  fluctuations  are  caused  random  impurity  scattering.  Shape- 
dependent  fluctuation  phenomena  are  a  fundamental  source  of  noise  disorder  for  any 
quantum  nanostructure  in  which  electrons  propagate  coherently. 

We  have  also  used  our  lithography  facilities  to  study  the  subband  structure  of 
electrons  in  wide  parabolic  quantum  wells  via  capacitance- voltage  (CV)  measurements. 
The  Fermi  energy  of  electrons  in  quantum  structures  introduces  a  quantum  term  to  the 
total  capacitance,  because  the  measured  voltage  between  an  electron  gas  and  a  gate  is 
actually  the  total  change  in  electrochemical  potential.  Rimberg  et  al.  [7]  demonstrate 
that  CV-profiling  measurements  on  parabolic  wells  with  in-plane  magnetic  fields  show 
features  which  can  be  unambiguously  associated  with  Fermi  level  changes  due  to  sub¬ 
band  depopulations,  in  addition  to  profiling  the  charge  density  of  the  electron  layer. 
In  the  future  we  plan  to  extend  the  sensitivity  of  this  technique  with  custom  made 
on  chip  transistors  to  do  quantum  CV  profiling  on  single  nanostructures,  as  recently 
demonstrated  by  Ashoori  and  coworkers. 

Nanostructure  electronic  devices  which  operate  via  quantum  mechanical  princi¬ 
ples  show  promise  for  applications  in  computing  and  signal  processing.  Although  their 
practical  use  lies  in  the  future,  we  can  study  the  physics  of  operation  of  nanostructure 
devices  now  in  order  to  determine  which  approaches  might  work  and  to  understand 
their  fundamental  limitations.  Because  the  temperature  of  operation  scales  inversely 
with  device  size,  what  we  learn  from  measurements  of  present  nanostructures  with 
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size  >  50  nm  at  very  low  temperatures  will  apply  to  smaller  structures  in  the  fu¬ 
ture  which  operate  at  more  accessible  temperatures.  We  have  excellent  facilities  for 
nanostructure  research  at  Harvard,  and  opportunities  for  novel  collaborative  efforts 
between  Prof.  Tinkham’s  group  (superconducting  nanostructures),  and  Prof.  Wester- 
velt’s  (semiconducting  nanostructures). 

We  propose  to  shift  our  efforts  toward  issues  associated  with  the  device  physics  of 
semiconductor  nanostructures  and  superconductor/semiconductor  hybrid  structures. 
Two  specific  projects  we  intend  to  pursue  are  the  development  of  custom-made  cryo¬ 
genic  field-effect  transistors  for  in  situ  measurements  of  single  nanostructures,  and  the 
development  of  a  Josephson  junction  oscillator  for  the  irradiation  of  semiconductor 
nanostructures  at  submillimeter  wavelengths. 

D  Cryogenic  Field  Effect  Transistors 

Ashoori  and  coworkers  have  recently  demonstrated  the  use  of  a  GaAs  heterostruc¬ 
ture  FET  to  make  capacitive  measurements  of  a  single  quantum  dot.  The  FET  is 
mounted  adjacent  to  the  dot  sample  and  both  sire  cooled  to  low  temperatures.  This 
configuration  is  naturally  suited  to  charge-coupled  measurements  of  nanostructures, 
because  the  input  resistance  is  very  high  and  the  input  capacitance  is  very  low  due 
to  the  small  size  of  the  gate  and  the  short  leads.  We  plan  to  extend  these  techniques 
by  constructing  matched  FET  pairs  using  the  two-dimensional  electron  gas  (2DEG) 
in  GaAs/AlGaAs  heterostructures.  We  currently  fabricate  gated  Hall  bar  structures 
routinely  for  quantum  Hall  effect  measurements,  using  either  optical  or  electron-beam 
lithography.  A  differential  FET  pair  can  be  constructed  from  two  matched,  gated 
2DEG  channels,  whose  input  capacitance  can  be  made  very  small  by  reducing  the 
size  of  the  transistors.  The  difference  in  charge  on  the  two  gates  can  be  measured 
with  high  accuracy  via  the  difference  in  channel  resistances  in  a  bridge  circuit.  The 
transistor  pair  can  be  fabricated  on  the  same  chip  as  the  nanostructure  sample  to 
achieve  extremely  low  input  capacitance,  or  on  a  separate  chip  for  measurements  in 
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a  magnetic  field  so  that  the  2DEG  channel  can  be  oriented  parallel  to  the  field  lines. 
The  extreme  sensitivity  possible  with  such  devices  should  permit  us  to  characterize 
single  nanostructures  in  a  wide  variety  of  experiments.  We  do  not  expect  difficulty  in 
producing  cryogenic  circuits,  because  the  yield  on  our  lithography  process  for  single 
devices  is  currently  quite  high. 

E.  Josephson  Junction/GaAs  Heterostructure  Hybrids 

Many  quantum  transitions  of  interest  for  semiconductor  heterostructures  and 
nanostructures  lie  in  the  submillimeter  wave  region  of  the  spectrum.  Spectroscopy  on 
nanostructures  is  difficult,  because  the  source  and  spectrometer  are  typically  located 
outside  the  sample  Dewar.  Channeling  submillimeter  wave  radiation  to  and  from  the 
sample  poses  special  problems,  especially  for  dilution  refrigerator  temperatures.  Fur¬ 
thermore,  arrays  are  usually  required,  because  single  nanostructures  typically  do  not 
provide  a  detectable  signal.  These  difficulties  can  be  avoided  by  constructing  a  super¬ 
conductor/semiconductor  hybrid  sample  in  which  a  Josephson  junction  is  used  as  a 
tunable  source  of  submillimeter  radiation.  By  placing  the  nanostructure  in  the  near 
field  of  the  oscillating  junction,  the  oscillating  electric  field  pattern  can  be  controlled, 
and  unwanted  radiation  into  the  far-field  can  be  largely  avoided,  so  that  all  leads  to  the 
sample  carry  only  dc  signals.  Hybrid  devices  such  as  these  are  naturally  suited  to  the 
measurement  of  single  nanostructures  at  dilution  refrigerator  temperatures.  In  collab¬ 
oration  with  Prof.  Tinkham’s  group,  we  have  already  begun  fabricating  test  structures 
in  which  a  Nb  Josephson  Junction  is  located  adjacent  to  a  quantum  point  contact  in 
a  GaAs  two-dimensional  electron  gas.  We  are  presently  working  to  improve  junction 
characteristics,  and  expect  results  in  the  near  future.  In  these  first  experiments  we 
intend  to  use  the  point  contact  as  a  threshold  detector  for  submillimeter  wave  pho¬ 
tons  to  demonstrate  the  viability  of  the  technique.  If  successful,  we  will  be  able  to  do 
spectroscopy  on  single  quantum  point  contacts  and  other  structures. 
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3.  R.  Seshadri  and  R.M.  Westervelt,  “Statistical  Mechanics  of  Magnetic  Bubble 
Arrays,  Part  1:  Topology  and  Thermalization,”  Phya.  Rev.  A  (to  appear). 

4.  R.  Seshadri  and  R.M.  Westervelt,  “Statistical  Mechanics  of  Magnetic  Bubble 
Arrays,  Part  2:  Observations  of  Two-Dimensional  Melting,”  Phya.  Rev.  A  (to 
appear). 

5.  R.  Seshadri  and  R.M.  Westervelt,  “Forced  Flow  of  Two-Dimensional  Magnetic 
Bubble  Arrays,”  Phya.  Rev.  Lett,  (to  be  submitted). 

6.  C.M.  Marcus,  A.J.  Rimberg,  R.M.  Westervelt,  P.F.  Hopkins,  and  A.C.  Gossard, 
“Conductance  Fluctuations  and  Chaotic  Scattering  in  Ballistic  Microstructures,” 
Phya.  Rev.  Lett,  (to  appear). 

7.  A.J.  Rimberg,  S.  Yang,  J.  Dempsey,  J.H.  Baskey,  R.M.  Westervelt,  M.  Sundaram, 
and  A.C.  Gossard,  “Capacitive  Depletion  of  the  Subband  Structure  in  the  Elec¬ 
tron  Gas  in  a  Wide  Parabolic  Quantum  Well,”  Appl.  Phya.  Lett,  (submitted). 
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ANNUAL  REPORT  OF 


PUBLICATIONS/PATENTS/PRESENTATIONS/HONORS 


a.  Papers  Submitted  to  Refereed  Journals  (and  not  yet  published) 

R.  Seshadri  and  R.M.  Westervelt,  ‘‘Two-Dimensional  Melting  of  Hexatic  Mag¬ 
netic  Bubble  Arrays,”  in  Proc.  Int.  Conf.  Two-Dimensional  Physics ,  Neuchatel, 
1991.  (Partially  supported  by  N00014-89-J-1592) 

R.  Seshadri  and  R.M.  Westervelt,  ‘‘Statistical  Mechanics  of  Magnetic  Bubble  Ar¬ 
rays,  Part  1:  Topology  and  Thermalization,”  Phys.  Rev.  A.  (Partially  supported 
by  N00014-89-J-1592) 

R.  Seshadri  and  R.M.  Westervelt.  "Statistical  Mechanics  of  Magnetic  Bubble 
Arrays,  Part  2:  Observations  of  Two-Dimensional  Melting,”  Phys.  Rev.  A.  (Par¬ 
tially  supported  by  N00014-S9-J-1592) 

R.  Seshadri  and  R.M.  Westervelt,  ‘‘Two- Dimensioned  Melting  of  Magnetic  Bubble 
Arrays:  A  Continuous  Hexatic  to  Liquid  Transition,”  Proc.  Fall  Meeting  of  the 
MRS,  Boston  1991.  (Partially  supported  by  N00014-89-J-1592) 

A.J.  Rimberg,  S.  Yang,  J.  Dempsey,  J.H.  Baskey,  R.M.  Westervelt,  M.  Sun- 
daram,  and  A.C.  Gossard,  “Capacitive  Depletion  of  the  Subband  Structure  in 
the  Electron  Gas  in  a  Wide  Parabolic  Quantum  Well,”  Appl.  Phys.  Lett.  (Par¬ 
tially  supported  by  N00014-89-J-1592) 

R.  Seshadri  and  R.M.  Westervelt,  “Forced  Flow  of  Two-Dimensional  Magnetic 
Bubble  Arrays,”  Phys.  Rev.  Lett.  (Partially  supported  by  N00014-89-J-1592) 

F.R.  Waugh  and  R.M.  Westervelt,  “Dynamics  of  Analog  Neural  Networks  with 
Local  Competition,”  Phys.  Rev.  A.  (Partially  supported  by  DARPA  AFOSR- 
89-0506) 

C.M.  Marcus,  A.J.  Rimberg,  R.M.  Westervelt,  P.F.  Hopkins,  and  A.C.  Gossard, 
“Conductance  Fluctuations  and  Chaotic  Scattering  in  Ballistic  Microstructures,” 
Phys.  Rev.  Lett.  (Partially  supported  by  NSF  DMR-88-17309  and  DARPA 
AFOSR-89-0506) 

b.  Papers  Published  in  Refereed  Journals 

C.M.  Marcus,  F.R.  Waugh,  and  R.M.  Westervelt,  “Connection  Topology  and 
Dynamics  in  Lateral  Inhibition  Networks,”  in  Advances  in  Neural  Information 
Processing  Systems  S  (Morgan  Kaufman,  San  Mateo,  1991),  p.  98.  (Partially 
supported  by  N00014-89-J-1592  and  AFOSR-89-0506) 
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C.M.  Marcias,  F.R.  Waugh,  and  R.M.  Westervelt,  “Nonlinear  Dynamics  and  Sta¬ 
bility  of  Analog  Neural  Networks,”  in  Nonlinear  Science :  The  Next  Decade ,  Proc. 
10th  Annual  Int.  Conf.  of  the  Center  for  Nonlinear  Studies,  Los  Alamos,  May 
1990,  D.  Campbell,  R.  Ecke,  and  J.M.  Hyman,  eds.  (North-Holland,  Amsterdam, 
1991),  p.  234.  (Partially  supported  by  N00014-89-J-1592) 

c.  Books  (and  sections  thereof)  Submitted  for  Publication 

F.R.  Waugh,  C.M.  Marcus,  and  R.M.  Westervelt,  “Nonlinear  Dynamics  of  Analog 
Associative  Memories,”  in  Neural  Associative  Memories,  ed.  M.  Hassoun  (Oxford 
University  Press,  Oxford,  1992).  (Partially  supported  by  DARPA  AFOSR-89- 
0506) 

g.  Invited  Presentations  at  Topical,  or  Scientific/Technical  Society  Con¬ 
ferences 

R.M.  Westervelt,  “Avalanches  in  Magnetic  Domain  Patterns,”  Santa  Fe  Institute, 
Workshop  on  Self- Organized  Criticality,  September  1991. 

R.M.  Westervelt,  “Space-Charge  Dynamics  in  Ultrapure  Ge,”  1st  Experimental 
Chaos  Conference,  Arlington,  Virginia,  October  1991.  (Partially  supported  by 
N00014-89-J-1592) 

R.M.  Westervelt,  “Two  Dimensional  Neural  Networks  for  Vision,”  Neural  Net¬ 
works  for  Computing,  Snowbird,  Utah,  April  1992. 

R.M.  Westervelt,  “Forced  Flow  of  Magnetic  Bubble  Arrays,”  12th  Annual  CNLS 
Conference  Nonlinearity  and  Materials  Modeling:  Simulations,  Interpretation 
and  Science,”  May  1992.  (Partially  supported  by  N00014-89-J-1592) 

R.M.  Westervelt,  “Chaotic  Scattering  of  Electrons  in  Quantum  Dots,”  Joint 
Soviet-American  Conference  on  Chaos,  Kiev,  Ukraine,  July  1992.  (Partially 
supported  by  N00014-89-J-1592) 

R.M.  Westervelt,  “Forced  Flow  of  Magnetic  Bubble  Arrays,”  Workshop  on  Spa¬ 
tially  Extended  Nonequilibrium  Systems,  Institute  for  Theoretical  Physics,  Uni¬ 
versity  of  California,  Santa  Barbara,  July  1992.  (Partially  supported  by  N00014- 
89-J-1592) 

C.M.  Marcus,  “Experimental  Investigations  of  Chaotic  Scattering  in  Semicon¬ 
ductor  Quantum  Dots,”  CHAOS  IV,  Kiev,  Ukraine,  July  1992. 

h.  Contributed  Presentations  at  Topical,  or  Scientific/Technical  Society 
Conferences 

C.M.  Marcus,  Meeting  of  the  American  Physical  Society,  Indianapolis,  Indiana, 
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March  1992. 


R.  Seshadri,  Meeting  of  the  American  Physical  Society,  Indianapolis,  Indiana, 
March  1992. 

R.  Seshadri,  International  Conference  on  Statistical  Physics,  STATPHYS  18, 
Berlin,  Germany,  August  1992. 

F.  Waugh,  International  Conference  on  Statistical  Physics,  STATPHYS  18, 
Berlin,  Germany,  August  1992. 

i.  Honors/Awards/Prizes 

R.M.  Westervelt,  Steering  Committee,  4th  Soviet- American  Chaos  Conference, 
1992. 

R.M.  Westervelt,  Naval  Studies  Board  Review  Panel,  Naval  Research  Laboratory, 
1992. 

R.M.  Westervelt,  Nonlinear  Dynamics  Program  Panel,  National  Science  Founda¬ 
tion,  1991. 

R.M.  Westervelt,  Member  of  JASON. 

j.  Graduate  Students  and  Postdoctoral  Supported  under  the  JSEP  for 
the  Year  Ending  July  31,  1902 


Dr.  Charles  Marcus,  Mr.  Fred  Waugh  and  Ms.  Raj  Seshadri 


1.7  Structural  and  Dynamical  Studies  of  Semiconductor  Interfaces  and 

Surfaces.  J.A.  Golovchenko,  Grant  N00014-89-J-1023;  Research  Unit  6. 

We  have  completed  the  structural  studies  of  lead  overlayers  on  silicon  and  ger¬ 
manium  surfaces  with  JSEP  support.  We  chose  these  two  systems  because  metal- 
semiconductor  surfaces  have  been  of  great  interest  in  science  as  well  as  in  technologi¬ 
cal  applications  for  decades.  Furthermore,  they  both  show  reversible  two-dimensional 
melting  transitions  which  merit  further  study.  Recently  it  was  found  on  several  systems 
(such  as  lead  on  silicon  and  nickel  on  silicon)  that  the  Schottky  barrier  height  depends 
strongly  on  the  interface  structure  upon  which  the  metal  overlayer  is  deposited.  Our 
studies  of  lead  on  silicon,  mainly  using  a  home-made  tunneling  microscope,  have  suc¬ 
cessfully  resolved  the  controversy  regarding  the  atomic  structures  and  the  number  of 
different  phases. 

Recently  we  extended  our  study  to  dynamical  processes  on  semiconductor  sur¬ 
faces.  Since  its  invention,  tunneling  microscopy  has  mainly  been  used  to  image  surface 
structures  at  atomic  resolution.  However,  its  use  in  the  study  of  dynamics  is  rare.  The 
study  of  atomic  scale  motions  is  of  fundamental  importance  in  the  understanding  of 
dynamical  phenomena  on  surfaces,  such  as  diffusion,  phase  transitions,  and  epitaxial 
growth.  Therefore,  our  tunneling  microscope  was  adapted  to  operate  at  temperatures 
above  25°C.  This  enables  us  to  observe  in  detail  motion  of  surface  atoms.  We  found 
that  a  small  number  of  lead  atoms  on  the  Ge(lll)  surface  can  catalyze  atomic  mo¬ 
tions  and  thus  allow  their  observation  near  room  temperature.  Our  initial  study  of  this 
system  showed  that  the  individual  interchanges  of  lead  and  germanium  atoms  within 
the  reconstructed  Ge(lll)-c(2  x  8)  surface  can  be  identified  and  that  the  interchange 
process  is  thermally  activated  and  obeys  the  Arrhenius  relation.  Its  activation  energy 
and  diffusion  constant  were  determined  by  a  plot  of  interchange  rate  vs.  1/T.  We 
also  found  that  this  diffusion  process  occurs  mainly  along  the  aidatom  row  direction  of 
the  c(2  x  8)  reconstructed  surface  and  that  about  hadf  of  the  interchanges  aire  “long 
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jumps”  (movements  of  more  than  one  atomic  spacing).  We  believe  this  is  the  first 
direct  observation  of  individual  diffusion  process  and  long  jumps.  A  paper  entitled 
“Direct  Measurement  of  Diffusion  by  Hot  Tunneling  Microscopy:  Activation  Energy, 
Anisotropy,  and  Long  Jumps”  was  published  in  March  1992  in  Physical  Review  Letters. 

Our  further  study  of  atomic  motions  for  lead  on  the  Ge(lll)  surface  showed  the 
formation  and  annihilation  of  metastable  structural  surface  excitations.  One  of  these 
excitations  is  associated  with  large  numbers  of  germanium  adatoms  shifting  positions 
along  the  adatom  row  like  beads  on  an  abacus.  This  effect  provides  a  new  mecha¬ 
nism  for  understanding  the  transport  of  atoms  on  a  semiconductor  surface.  Also,  it 
earn  account  for  the  previously  observed  anisotropic  diffusion  and  provides  clues  for 
the  phase  transition  of  pure  Ge(lll)-c(2  x  8)  to  1  x  1  phase,  which  is  still  not  well 
understood.  We  believe  our  observation  of  atomic  motions  will  stimulate  further  the¬ 
oretical  calculations  which  may  enhance  the  fundamental  understanding  of  dynamical 
processes.  A  manuscript  entitled  “Observation  of  Metastable  Structural  Surface  Ex¬ 
citations  and  Concerted  Adatom  Motions:  Pb/Ge(lll)”  was  submitted  recently  to 
Science  magazine.  A  full  report  of  our  observations  of  atomic  motions  on  this  surface 
is  in  preparation. 

A  second  tunneling  microscope  has  been  constructed  in  our  laboratory.  Clear 
images  of  Si(lll)-7  x  7  have  been  obtained  at  room  temperature;  however,  minor  mod¬ 
ifications  of  the  microscope  and  of  the  control  electronics  and  software  are  still  needed. 
This  STM  is  also  designed  so  that  samples  may  be  held  at  elevated  temperature  during 
tunneling,  to  facilitate  additional  studies  of  diffusion  and  growth  kinetics.  It  will  be 
used  to  continue  our  studies  of  diffusion  on  semiconductor  surfaces  and  to  study  the 
growth  of  Ge  on  Si(lll). 

We  have  initiated  studies  of  the  structural  and  electronic  properties  of  hydrogen- 
and  halogen- terminated  silicon  (111)  surfaces.  It  is  known  that  atomically  flat  and 
electronically  passivated  hydrogen- terminated  silicon  (111)  surfaces  can  be  obtained 
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easily  from  chemical  etching  of  a  silicon  wafer  in  buffered  hydrofluoride  solution.  Each 
surface  silicon  atom  is  bound  by  one  hydrogen  atom,  which  saturates  the  surface 
dangling  bond  and  passivates  the  surface.  On  the  other  hand,  a  reliable  procedure  to 
produce  ideal  halogen- terminated  silicon  surfaces  is  still  to  be  found.  To  determine 
the  surface  quality  of  trial  samples  of  halogen- terminated  silicon  surfaces,  an  X-ray 
standing  wave  setup  has  been  completed  and  is  being  used  to  measure  the  position 
of  bromine  atoms  chemisorbed  on  silicon  (111)  surfaces.  Preliminary  experimental 
results  are  consistent  with  the  picture  that  one  bromine  atom  is  bound  to  one  surface 
silicon  atom.  The  measured  distance  between  bromine  atoms  and  bulk  silicon  (111) 
plane  is  2.57  A. 

The  mechanism  of  surface  passivation  and  the  formation  of  Schottky  barriers  are 
two  important  subjects  in  solid  state  electronics.  We  will  extend  our  study  of  hydrogen- 
and  halogen-terminated  silicon  surfaces  to  these  two  subjects  and  look  for  correlations 
between  the  structured  and  electronic  properties  of  the  surfaces.  Stimulated  by  the 
realization  of  organic  electroluminescent  devices,  we  are  looking  into  the  possibility 
of  binding  dye  molecules  on  silicon  surfaces  and  making  light-emitting  devices  out  of 
them. 
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ANNUAL  REPORT  OF 

PUBLICATIONS /PATENTS /PRESENTATIONS/HONORS 


a.  Papers  Submitted  to  Refereed  Journals  (and  not  yet  published) 

I. -S.  Hwang  and  J.A.  Golovchenko,  ‘‘Observation  of  Metastable  Structural  Sur¬ 
face  Excitations  and  Concerted  Adatom  Motions:  Pb/Ge(lll),”  Science. 

b.  Papers  Published  in  Refereed  Journals 

E.  Ganz,  F.  Xiong,  I.-S.  Hwang,  and  J.A.  Golovchenko,  “Submonolayer  Phases 
of  Pb  on  Si(lll),”  Phys.  Rev.  B43,  7316  (1991). 

E.  Ganz,  I.-S.  Hwang,  F.  Xiong,  S.K.  Theiss,  and  J.A.  Golovchenko,  “Growth 
and  Morphology  of  Pb  on  Si(lll),”  Surface  Science  257,  259  (1991). 

E.  Ganz,  S.K.  Theiss,  I.-S.  Hwang,  and  J.A.  Golovchenko,  “Direct  Measurement 
of  Diffusion  by  Hot  Tunneling  Microscopy:  Activation  Energy,  Anisotropy,  and 
Long  Jumps,”  Phys.  Rev.  Lett.  68,  1567  (1992). 

g.  Invited  Presentations  at  Topical,  or  Scientific/Technical  Society  Con¬ 
ferences 

J. A.  Golovchenko,  “Metastable  Structural  Surface  Excitations  and  Correlated 
Atomic  Motion:  A  STM  Study  of  the  Mechanism  of  Surface  Diffusion  within  a 
Semiconductor  Surface,”  March  1992  Meeting  of  the  American  Physical  Society. 

h.  Contributed  Presentations  at  Topical,  or  Scientific/Technical  Society 
Conferences 

E.  Ganz,  F.  Xiong,  I.-S.  Hwang,  and  J.A.  Golovchenko,  “Growth  and  Morphology 
of  Pb  on  Si(lll)  7  x  7,”  Abstract  for  the  18th  Annual  Conf.  on  the  Physics  and 
Chemistry  of  Semiconductor  Interfaces ,  Long  Beach,  California  (1991). 

j.  Graduate  Students  and  Postdoctorals  Supported  under  the  JSEP  for 
the  Year  Ending  July  SI,  1992 

Ms.  Silva  K.  Theiss  and  Mr.  Ing-Shouh  Hwang 
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II.  QUANTUM  ELECTRONICS 


Personnel 


Prof.  E.  Mazur 
Dr.  W.  Mieher 
Mr.  A.  Feder 
Mr.  J.  Goldman 
Mr.  S.  Deliwala 


Ms.  K.Y.  Lee 
Mr.  C.Z.  Lu 
Mr.  E.  Glezer 
Mr.  Y.  Siegal 
Mr.  J.K.  Wang 


H.l  Femtosecond  Nonlinear  Optical  Interactions  E.  Mazur,  E.  Glezer,  Y.  Sie¬ 
gal,  and  J.K.  Wang,  Grants  N00014-89-J-1023  and  NSF  DMR-8858075;  Research 
Unit  7. 

In  this  past  year  we  have  taken  extensive  measurements  of  the  response  of  GaAs 
to  excitation  with  high-energy,  150-fs  laser  pulses.  When  excited  above  a  threshold 
fluence  of  100  mJ/cm2,  the  electronic  system  in  the  material  undergoes  an  ultrafast 
transformation  to  a  disordered,  metallic  phase.  The  disordering  is  evidenced  by  the 
vanishing  of  the  reflection  second  harmonic  signal  of  a  probe  pulse  with  a  100-fs  decay 
time  while  a  sharp  rise  in  linear  reflectivity  indicates  the  metallic  character  of  the 
disordered  phase.  These  results  suggest  that  lattice  disordering  can  be  driven  directly 
by  the  breaking  of  bonds  through  laser  excitation  of  valence  electrons  to  the  conduction 
band  and  that  this  disordering  can  occur  while  the  lattice  is  still  cold. 

In  addition  to  the  disordering  measurements,  we  have  obtained  preliminary  re¬ 
sults  on  ultrahigh  density  carrier  recombination  below  the  disordering  threshold  using 
a  highly  sensitive  probe  of  the  index  of  refraction  involving  reflection  second  harmonic 
generation.  Little  is  known  about  the  behavior  of  GaAs  in  the  high  carrier  density 
regime,  and  our  results  show  an  Auger  carrier  recombination  time  that  is  significantly 
faster  than  that  predicted  by  theory.  At  carrier  densities  near  1022/cm3,  the  measured 
recombination  time  is  about  500  fs  instead  of  the  10  ps  value  predicted  by  a  static 
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screening  model. 

We  recently  completed  rebuilding  the  laser  amplifier  with  a  novel,  flexible  design 
involving  continuum  generation  in  a  single-mode  optical  fiber  followed  by  amplifica¬ 
tion  in  prism  dye  cells.  This  design  allows  us  to  have  two  high-energy,  independently 
tunable,  synchronized  beams  of  60-fs  laser  pulses  with  excellent  spatial  profile.  With 
this  added  flexibility  we  are  now  in  a  position  to  embark  on  a  vigorous  set  of  experi¬ 
ments  designed  to  fully  explore  ultrafast  laser-induced  disordering  as  well  as  ultrahigh 
density  carrier  dynamics  in  semiconductors. 

(i)  Studies  of  carrier  dynamics  at  very  high  densities.  Using  our  sensitive  re¬ 
flection  second  harmonic  probe  of  the  index  of  refraction,  we  plan  to  examine  more 
closely  the  behavior  of  carriers  in  crystalline  GaAs  at  conduction  electron  densities 
greater  than  1021.  We  can  easily  achieve  these  carrier  densities  with  our  femtosecond 
laser  pulses  because  the  laser  energy  is  deposited  on  such  a  short  time  scale.  Very 
little  information  exists  on  this  regime  of  carrier  dynamics,  yet  it  is  becoming  more 
and  more  important  to  understand  the  physics  at  high  carrier  densities  because  this  is 
precisely  the  situation  in  many  new  high  technology  devices.  For  example,  because  of 
the  gain  limiting  effects  of  recombination,  information  on  high  carrier  density  Auger 
recombination  times  will  be  crucial  in  the  rapidly  developing  field  of  semiconductor 
diode  lasers. 

(ii)  Temperature  dependence  of  semiconductor  disordering.  Since  the 
electronically-induced  disordering  occurs  before  significant  energy  is  transferred  from 
the  hot  electrons  to  the  lattice,  the  temperature  of  the  lattice  upon  disordering  is 
probably  not  much  higher  than  room  temperature.  The  average  thermal  velocity  of 
the  ions  at  room  temperature  is  fast  enough  to  lead  to  a  metallic,  disordered  lattice 
in  a  few  hundred  femtoseconds  if  the  bonds  that  hold  the  ions  in  place  are  suddenly 
removed.  If  this  bond-breaking  picture  is  accurate,  then  the  measured  disordering 
time  should  depend  on  the  average  thermal  velocity  of  the  ions  upon  excitation  with 
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the  laser  pulse.  We  plan  to  build  a  cryogenic  system  that  will  allow  us  to  cool  the 
sample  down,  to  below  20°  K.  This  will  allow  us  to  study  the  effects  of  slower  initial  ion 
velocities  on  the  disordering  time. 

(iii)  Ellipsometry  of  Disordered,  Semiconductors.  Although  we  know  that  the 
short-lived  phase  of  GaAs  resulting  from  electronically-induced  disordering  is  metallic 
in  character,  many  more  details  remain  to  be  learned  about  its  nature.  The  dielectric 
constant  at  620  nm  which  we  extracted  from  our  data  does  not  match  with  theoretical 
simulations  for  the  dielectric  constant  of  equilibrium  liquid  GaAs,  suggesting  that  fem¬ 
tosecond  laser-induced  disordering  may  result  in  a  state  unattainable  through  other 
means.  By  mapping  out  the  index  of  refraction  of  disordered  semiconductors  using 
pump-probe  ellipsometry,  we  can  obtain  a  fairly  detailed  picture  of  these  novel  mate¬ 
rials.  Our  new,  tunable,  laser  amplifier  gives  us  the  ability  to  study  the  wavelength 
dependence  of  the  index  of  refraction. 

n.3  Interactions  with  Ultrashort  Laser  Pulses.  E.  Mazur,  C.Z.  Lu,  S.  Deli- 

wala,  and  J.  Goldman,  Grants  N00014-89-J-1023  and  ARO  DAAL03-88-K-0114; 

Research  Unit  8. 

In  this  research  unit  we  made  significant  progress  this  year  towards  our  goal 
of  developing  an  experimental  apparatus  for  femtosecond  laser  studies  of  molecular 
surface  processes  in  ultrahigh  vacuum  (UHV).  We  completed  renovation  of  our  clean 
air  lab  facility,  constructed  a  femtosecond  Ti:Sapphire  laser  oscillator,  designed  a  high- 
power  high-repetition-rate  laser  amplifier,  and  designed  and  purchased  a  UHV  surface 
science  chamber. 

The  self-mode-locked  TitSapphire  laser  oscillator  is  currently  operational  and 
generates  pulses  of  80-fs  duration.  We  chose  Ti:Sapphire  because  of  its  great  versatility 
as  a  short  pulse  source,  being  tunable  with  a  single  set  of  optics  from  700-900  nm.  The 
laser  has  excellent  mechanical  stability.  A  HITCI  saturable  absorber  jet  in  the  cavity 
stabilizes  the  self-mode  locking  process  and  provides  long  term  pulse-train  stability 
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required  for  our  high-repetition  rate  multi-shot  averaging  techniques.  A  surprising 
advantage  of  the  Ti:Sapphire  oscillator  is  that  it  can  produce  very  high  average  powers. 
Using  a  20%  output  coupler,  we  obtain  several  watts  of  mode-locked  infrared  output 
for  10-15  W  of  argon  laser  pump  power.  This  means  that  we  can  obtain  tens  of  nJ  per 
laser  pulse,  which  is  sufficient  for  continuum  generation  or  second  harmonic  probing 
of  surfaces  without  further  amplification.  These  high  pulse  energies  also  eliminate  the 
need  for  multiple  stages  of  amplification;  with  a  simple  six-pass  amplifier  design  we  will 
be  able  to  obtain  about  1-10  fiJ  per  pulse,  and  additional  linear  amplifier  stages  will 
allow  us  to  produce  energies  in  the  100  pJ  range  if  we  desire  such  pulses  for  disordering, 
desorption  or  reaction  studies.  Currently  we  are  designing  the  amplifier  section  of  the 
laser.  Until  funds  can  be  obtained  for  a  high  repetition  rate  (>  1  kHz)  pump  laser, 
we  plan  to  use  our  existing  10-Hz  Nd:YAG  laser  in  a  linear  amplifier  design.  This  will 
allow  us  to  obtain  preliminary  data;  on  the  long-term  a  high-repetition  rate  laser  with 
a  regenerative  amplifier  would  allow  us  to  increase  the  signal- to-noise  ratio  by  several 
orders  of  magnitude. 

Our  proposed  laser  surface  studies  in  UHV  require  a  mechanically  stable,  com¬ 
pact  chamber  with  short  optical  path  length  access  to  the  sample  over  a  wide  range  of 
angles.  To  this  end  we  designed  a  two  level  chamber,  with  a  vertically  mounted  VG 
manipulator  which  allows  precision  sample  positioning  in  the  upper  or  lower  cham¬ 
bers.  The  lower  level  is  a  small  six  inch  cube  with  windows  on  four  sides,  allowing 
positioning  of  lenses  or  other  optical  components  outside  the  chamber  only  two  inches 
away  from  the  sample;  with  the  sample  in  this  position  we  can  achieve  a  continuous  45 
degree  field  of  view  which  makes  the  alignment  of  several  beams  in  a  phase  matched 
geometry  much  more  practical.  Backing  the  sample  up  to  four  inches  from  the  windows 
allows  grazing  incidence  geometries  as  well.  The  upper  chamber  level  is  equipped  with 
standard  surface  preparation  tools,  a  load  lock,  and  characterization  probes  including 
mass  spectrometry  and  LEED/Auger  spectroscopy.  The  manipulator  is  designed  to 
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allow  rapid  sample  cooling  to  IN?  temperatures  (down  to  100  K)  which  is  required  for 
preparation  of  oxygenated  metal  surfaces. 

The  most  immediate  goal  for  this  apparatus  is  to  time-resolve  a  two-species  laser- 
induced  surface  chemical  process.  In  particular  we  are  interested  in  the  properties  of 
metal  catalytic  surfaces  with  monolayer  adsorbed  coverages  of  gaseous  molecules  such 
as  oxygen.  We  hope  to  test  severed  nonlinear  optical  techniques  for  their  effectiveness 
in  characterizing  these  surfaces  as  they  undergo  molecular  steps  of  reaction  including 
vibrational  and  electronic  excitation  and  relaxation,  dissociation  and  desorption,  and 
recombination  into  new  species.  Such  studies  offer  enormous  potential  for  furthering 
the  understanding  of  these  fundamental  molecular  processes,  and  for  the  possibility  of 
laser-driven  final-state-specific  laser  chemistry  and  materials  processing.  Experiments 
leading  up  to  these  goals  will  include  laser-induced  dissociation  of  O2  and  detection 
of  OH  from  a  laser  excited  O2/H2  coadsorbate  surface.  In  addition  to  the  surface 
reaction  work,  we  plan  to  investigate  the  effect  of  laser-excited  adsorbates  on  the 
electrical  properties  of  semiconductor  and  met  ad  surfaces  —  a  study  that  may  have 
significant  import  for  device  physics  applications. 
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J. -K.  Wang,  Y.  Siegal,  C.  Lii,  and  E.  Mazur,  “‘Generation  of  Synchronized,  Tun¬ 
able,  Femtosecond  Laser  Pulses  with  Nearly  Perfect  Gaussian  Spatial  Profile,” 
Opt.  Comm. 

K. -H.  Chen,  C.-Z.  Lii,  J.  Goldman,  S.  Deliwala,  and  E.  Mazur,  “Coherent  Anti- 
Stokes  Raman  Spectroscopy  of  Infrared  Multiphoton  Excited  Molecules,”  J. 
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P.  Saeta,  Y.  Siegal,  J.-K.  Wang,  and  E.  Mazur,  “Femtosecond  Carrier-induced 
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P.  Saeta,  J.-K.  Wang,  Y.  Siegal,  N.  Bloembergen,  and  E.  Mazur,  “Ultrafast 
Electronic  Disordering  During  Femtosecond  Laser  Melting  of  GaAs,”  Phys.  Rev. 
Lett.  67,  1023  (1991). 

C.-Z.  Lii,  J.  Goldman,  S.  Deliwala,  K.-H.  Chen,  and  E.  Mazur,  “Direct  Evidence 
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J.-K.  Wang,  Y.  Siegal,  P.  Saeta,  N.  Bloembergen,  and  E.  Mazur,  “Femtosec¬ 
ond,  Electronically- Induced  Disordering  oi  GaAs,”  in  Ultrafast  Phenomena  VIII, 
(Springer  Verlag,  1992)  (in  press). 

E.  Mazur,  J.-K.  Wang,  and  Y.  Siegal,  “Femtosecond  Nonlinear  Optics  of  Semi¬ 
conductors,”  Proe.  Army  Sei.  Conf.  (in  press). 

d.  Books  (and  sections  thereof)  Published 

P.  Saeta,  J.-K.  Wang,  Y.  Siegal,  N.  Bloembergen,  and  E.  Mazur,  “Electronic 
Disordering  During  Femtosecond  Laser  Melting  of  GaAs,”  in  Technical  Digest, 
Quantum  Electronics  and  Laser  Science  Conference,  Anaheim,  CA,  May  1992. 
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K.Y.  Lee,  T.  Chou,  and  E.  Mazur,  “Direct  Optical  Measurements  of  Capillary 
Wave  Damping  at  Liquid- Vapor  Interfaces,”  Technical  Digest,  Int.  Quantum 
Electronics  Conf.,  Vienna,  June  1992. 

J.-K.  Wang,  Y.  Siegal,  C.  Lii,  and  E.  Mazur,  “Subpicosecond  Stimulated  Raman 
Scattering  in  High-pressure  Hydrogen,”  Technical  Digest,  Int.  Quantum  Elec¬ 
tronics  Conf.,  Vienna,  June  1992. 

g.  Invited  Presentations  at  Topical  or  Scientific /Technical  Society  Con¬ 
ferences 

E.  Mazur,  “Electronic  Disordering  During  Femtosecond  Melting  of  Gallium  Ar¬ 
senide,”  Quantum  Electronics  and  Laser  Symposium,  Anaheim,  CA,  May  1992. 

E.  Mazur,  Keynote  Address:  “Femtosecond  Nonlinear  Optics  of  Semiconduc¬ 
tors,”  Army  Science  Conference,  Orlando,  Florida,  June  1992. 

h.  Contributed  Presentations  at  Topical  or  Scientific/Technical  Society 
Conferences 

E.  Mazur,  “Reflectivity  and  Second  Harmonic  Efficiency  of  GaAs  During  Fem¬ 
tosecond  Melting,”  Ultrafast  Phenomena  Conference,  Monterey,  California,  June 
1990. 

• 

E.  Mazur,  “Collisional  and  Intramolecular  Dynamics  of  Low  Lying  Vibrational 
States  of  Infrared  Multiphoton  Excited  Molecules,”  International  Quantum  Elec¬ 
tronics  Conference,  Anaheim,  California,  June  1990. 

E.  Mazur,  “Dynamics  of  Vibrational  Energy  Transfer  in  IR  Laser-excited 
Molecules  Observed  Using  Nonlinear  Spectroscopy,”  Interdisciplinary  Laser  Sci¬ 
ence  Conference  1990,  Minneapolis,  September  1990. 

E.  Mazur,  “Femtosecond  Melting  of  GaAs  Probed  by  Reflectivity  and  Second- 
harmonic  Generation,”  1990  Annual  Meeting  of  the  Optical  Society  of  America, 
Boston,  November  1990. 

E.  Mazur,  “Tunable,  Synchronized,  High  Energy,  Femtosecond  Laser  Pulses  with 
Gaussian  Spatial  Profile,”  1991  Annual  Meeting  of  the  Optical  Society  of  Amer¬ 
ica,  San  Jose,  November  1991. 

E.  Mazur,  “Femtosecond,  Electronically-Induced  Disordering  of  GaAs,”  Ultrafast 
Phenomena  Conference,  Antibes,  France,  June  1992. 

E.  Mazur,  “Subpicosecond  Stimulated  Reiman  Scattering  in  High-pressure  Hy¬ 
drogen,”  International  Conference  on  Quemtum  Electronics  ’92,  Vienna,  Austria, 
June  1992. 
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E.  Mazur,  “Direct  Optical  Measurements  of  Capillary  Wave  Damping  at  Liquid- 
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CIRCUITS  AND  SYSTEMS 


Personnel 


Prof.  R.W.  Brockett 

Assoc.  Prof.  J.J.  Clark 

Prof.  Y.C.  Ho 

Dr.  L.  Mollamustafaoglu 

Dr.  R.  Sreenivas 

Dr.  Z.  B.  Tang 

Asst.  Prof.  W.  Yang 

Mr.  C.H.  Chen 

Mr.  L.-Y.  Dai 


Mr.  M.  Deng 
Mr.  D.  Friedman 
Mr.  I.  Garai 
Mr.  M.  Gottardi 
Mr.  W.G.  Li 
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Mr.  N.  Saxena 
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m.l  Analog  VLSI  and  Principal  Component  Analysis.  R.W.  Brockett, 

J.J.  Clark,  W.  Yang,  and  N.  Saxena,  Grant  N00014-89-J-1023;  Research 

Unit  9. 

Since  November  1991,  we  have  been  investigating  adaptive,  nonlinear,  signal 
processing  tasks  and  have  focussed  on  the  task  of  principal  component  analysis. 
Principal  component  analysis  is  a  basic  signal  processing  task  that  has  important 
applications  in  data  compression,  noise  rejection,  and  pattern  recognition.  In  recent 
work  by  Brockett  at  Harvard  University,  we  have  described  a  system  of  nonlinear 
differential  equations  which  can  determine  the  eigenvalues  of  a  sample  covariance 
and  thus  be  used  for  principal  component  analysis.  While  previous  VLSI  imple¬ 
mentations  of  these  types  of  adaptive  signal  processing  tasks  have  utilized  digital 
components,  the  characteristics  of  coupled  differential  equations  and  principal  com¬ 
ponent  analysis  are  well  matched  to  the  capabilities  of  analog  VLSI.  Currently, 
we  are  investigating  alternative  VLSI  implementations  which  utilize  novel  analog 
circuitry  and  massively  parallel  architectures.  A  robust  analog  computing  system 
offers  potentially  revolutionary  improvements  in  speed,  size,  power  consumption, 
and  cost  in  comparison  to  traditional  digital  approaches.  However,  analog  VLSI 
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implementations  must  be  able  to  deal  effectively  with  noise,  limited  dynamic  range, 
imprecise  component  values,  mismatched  components,  and  accumulation  of  errors 
both  at  the  circuitry  and  algorithmic  level.  Our  research  efforts  have  focussed  on 
addressing  these  issues  for  the  implementation  of  a  robust  analog  computing  system 
for  principal  component  analysis. 

Specifically,  principal  component  analysis  involves  the  computation  of  an 
N  x  N  sample  covariance  of  a  signal  and  subsequently  determining  the  associated 
N  eigenvalues  of  the  sample  covariance.  Through  a  series  of  multiply  and  accu¬ 
mulate  operations,  we  can  compute  the  sample  covariance  matrix.  This  sample 
covariance  matrix  is  used  as  the  initial  conditions  for  our  system  of  nonlinear  differ¬ 
ential  equations  which  will  evolve  in  time  to  yield  the  eigenvalues  which  are  sorted 
from  largest  to  smallest.  The  corresponding  eigenvectors  can  also  be  found  using  a 
related  system  of  equations.  Thus,  a  full  implementation  of  a  principal  component 
analyzer  will  require  two  components,  one  for  computing  the  sample  covariance  and 
one  for  implementing  our  system  of  nonlinear  differential  equations. 

III.2  Analog  VLSI  Implementation  of  the  Eigenvalue  Finder.  R.W. 
Brockett,  J.J.  Clark,  W.  Yang,  and  N.  Saxena,  Grant  N00014-89-J-1023; 

Research  Unit  9. 

An  important  aspect  of  the  implementation  of  the  dynamical  systems  which 
perform  the  matrix  diagonalization  is  that  of  the  effect  of  circuit  inaccuracies  on 
the  convergence  of  the  system.  Examination  of  the  equations  indicates  that,  de¬ 
pending  on  the  particular  approach  taken,  our  system  requires  the  implementation 
of  analog  multipliers  or  exponentiation  circuits,  scaling  circuits  (i.e.,  multiplication 
by  a  constant),  adders  and  integrators. 

All  these  component  circuits  are  quite  common  in  the  analog  computing 
world.  However,  as  is  quite  well-known,  such  circuits  have  limited  accuracy  (about 
6  to  8  bits  at  best)  which  could  severely  affect  the  performance  of  the  entire  com- 
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putation.  It  is  crucial  that  our  implementation  of  the  system  of  nonlinear  equations 
behaves  properly  in  the  face  of  this  level  of  precision.  Another  complication  is  the 
limited  dynamic  range  of  analog  computing  circuits,  which  could  also  affect  the  per¬ 
formance  of  the  computing  algorithm.  In  order  to  get  an  idea  of  the  effect  of  limited 
accuracy  and  dynamic  range  on  the  entire  computation,  a  simulation  experiment 
was  performed.  A  system  of  3  x  3  equations  was  simulated  using  MATLAB.  Lim¬ 
ited  accuracy  was  introduced  by  simply  truncating  the  result  of  each  computation 
after  the  second  decimal  bit.  Limited  dynamic  range  was  introduced  by  keeping  the 
result  of  each  computation  within  a  certain  range  of  values  (for  example,  between 
0.01  and  3  volts).  The  input  values  were  restricted  between  0.01  and  1  volts.  The 
result  obtained  from  simulations  with  these  inaccuracies  was  then  compared  with 
that  obtained  from  performing  the  same  computation  using  32  bit  floating  point 
operations.  The  performance  of  the  system  with  the  inaccuracies  was  qualitatively 
similar  to  that  of  the  system  without  the  inaccuracies. 

Two  approaches  to  the  implementation  of  the  equations  are  being  exam¬ 
ined.  The  first  involves  the  use  of  analog  multipliers  while  the  second  implements 
a  transformed  version  of  the  equations  which  require  exponentiators  and  summa¬ 
tion  circuits  in  place  of  the  multipliers.  Currently,  we  are  working  on  the  second 
approach  and  will  begin  work  on  the  first  approach  in  the  near  future. 

The  approach  using  exponentiators  is  promising  as  the  operation  of  expo¬ 
nentiation  (and  consequently  of  logarithms)  can  be  done  quite  naturally  and  simply 
in  VLSI  technology.  Two  approaches  to  its  implementation  are  being  considered: 
i)  MOS  Field  Effect  Transistors  (this  approach  uses  the  exponential  characteristic 
of  a  MOSFET  in  its  subthreshold  regime),  ii)  Bipolar  junction  (PN  Junction  Diode 
and  Parasitic  Bipolar  Transistor). 

A  test  chip  with  several  bipolar  transistor  configurations  and  an  MOS  expo¬ 
nential  circuit  has  been  fabricated  and  is  currently  being  tested.  Once  we  know  the 
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exponential  characteristics  of  the  fabricated  devices  and  circuits,  we  will  be  able  to 
design  the  multiplier  more  accurately. 

Simulations  show  that  the  bipolar  junction  has  a  larger  range  of  exponential 
behavior  than  the  MOS  circuit.  So,  while  the  chip  is  being  tested,  the  bipolar 
junction  transistor  is  being  used  for  design  simulations  of  the  multiplier  circuit. 
The  logarithm  circuit  was  designed  using  an  opamp  and  a  bipolar  transistor.  The 
output  of  two  such  circuits  was  fed  into  the  inputs  of  an  opamp  adder,  whose  output 
was  then  exponentiated  using  a  single  bipolar  transistor.  Thus  we  get  a  current 
output  which  is  a  scaled  multiplication  of  the  current  inputs.  This  design  was 
simulated  and  was  found  to  work  for  a  suitable  range  of  input  values.  The  results 
have  not  been  analyzed  in  much  detail  because  the  model  for  the  bipolar  transistor 
used  in  the  simulation  is  not  for  parasitic  bipolar  transistors.  Once  the  parameters 
for  the  parasitic  bipolar  transistor  are  obtained,  a  more  realistic  simulation  will 
be  performed,  and  the  simulation  results  will  be  analyzed  in  more  detail.  One  can 
avoid  the  need  for  logarithm  operations  if  the  logarithms  of  the  entries  of  the  matrix 
to  be  diagonalized  are  computed  externally.  Then  the  matrix  differential  equations 
need  only  require  exponentiation  operations. 

III.3  Toda  Lattice  Sorter.  R.W.  Brockett,  J.J.  Clark,  W.  Yang,  and  N.  Sax- 

ena,  Grant  N00014-89-J-1023;  Research  Unit  9. 

An  interesting  special  case  of  the  symmetric  matrix  diagonalization  process 
occurs  when  the  matrix  to  be  diagonalized  is  diagonal  but  with  its  diagonal  entries 
not  sorted.  The  matrix  differential  equations  will  flow  in  such  a  way  as  to  move  away 
from  the  initial  condition  and  then  return  to  a  diagonal  matrix  in  which  the  diagonal 
entries  are  sorted.  It  can  be  shown  that  the  matrix  is  at  all  times  tri-diagonal,  which 
greatly  simplifies  the  connectivity  required  for  its  implementation.  In  fact  one  can 
implement  this  sorter  (known  as  a  “Toda”  lattice)  with  a  one- dimensional  array  of 
processing  elements. 
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We  have  fabricated  a  CMOS  integrated  circuit  which  contains  a  30-stage 
Toda  lattice  sorter  as  well  as  some  test  circuitry.  This  circuit  uses  the  exponential 
form  of  the  equations  and  will  give  us  data  on  the  performance  of  the  approaches  we 
are  taking  to  the  implementation  of  exponentiation.  We  have  very  recently  received 
the  fabricated  circuits  and  will  begin  characterizing  their  performance  in  the  near 
future. 


m.4  Massively  Parallel  Processor  for  Eigenvalue  Computation.  R.W. 

Brockett,  J.J.  Clark,  W.  Yang,  and  N.  Saxena,  Grant  N000 14-89- J-l 023; 

Research  Unit  9. 

As  has  been  previously  shown,  the  eigenvalues  of  a  symmetric,  positive  def¬ 
inite  matrix,  H( 0)  can  be  found  as  the  diagonal  elements  of  the  solution,  H( oo) 
to 

H  =  [H,[H,N\] 

where  [A,  B]  —  AB  -  BA  and  N  is  chosen  to  be  diagonal  with  nlt  =  i.  Thus  for 
each  element  in  the  matrix,  h,}  we  can  write 

n 

htj  =  ^  ^  hikhkj{i  "b  j  Ik). 

k=  1 

Notice  that  the  computation  of  ht]  requires  the  weighted  sum  of  a  series  of  analog- 
analog  multiplications.  Furthermore,  notice  that  the  computation  of  hij  does  not 
require  arbitrary  communication  capability  but  instead  only  requires  communica¬ 
tion  within  common  column  elements,  /i,*  and  row  elements,  hkj. 

By  exploiting  this  limited  communications  requirement,  we  can  implement 
a  massively  parallel  processor  for  performing  principal  component  extraction  (see 
Figure  III.  1 ).  The  elements,  htJ  of  a  symmetric,  positive  definite  matrix  are  loaded 
into  the  processor  array  as  indicated.  Then  each  value,  htJ  is  broadcast  to  common 
column  and  row  elements  through  vertical  and  horizontal  connections.  The  products 
of  column  and  row  elements  axe  computed  by  processors  at  each  location  indicated 
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Figure  lll.l.  A  massively  parallel  processor  for  computing  eigenvalues  of  a  4  x  4  matrix. 

by  a  circle.  Thus,  the  value  of  hij  is  computed  as  the  indicated  weighted  sum  of 
the  products  of  row  and  column  elements.  When  hX)  =0  for  every  element  in  the 
array,  the  system  has  reached  steady  state,  the  off-diagonal  elements  will  be  0,  and 
the  eigenvalues  will  be  ordered  from  largest  to  smallest  as  the  diagonal  elements, 
h{t. 


in.5  Pulse-Frequency  Analog  Computing.  R.W.  Brockett,  J.J.  Clark,  W. 
Yang,  D.J.  Friedman,  and  N.  Saxena,  Grant  N00014-89-J-1023;  Research 
Unit  9. 

Analog  information  has  typically  been  encoded  as  an  analog  change,  cur¬ 
rent,  or  voltage  value  which  effectively  limits  the  dynamic  range  in  simple  im¬ 
plementations  to  103  due  to  noise  and  signail  degradation.  Recently,  we  developed 
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and  demonstrated  a  pulse-frequency  technique  for  robust  encoding  of  wide-dynamic 
range  analog  information.  In  particular,  we  implemented  a  VLSI  array  of  32  x  32 
photosensors  with  a  linear  dynamic  range  of  107  in  optical  energy.  By  encoding 
analog  information  in  the  time  domain  or  frequency  domain,  pulse-frequency  en¬ 
coding  is  much  less  susceptible  to  degradation  and  can  extend  the  dynamic  range  of 
analog  VLSI  components.  We  are  currently  utilizing  this  pulse-frequency  technique 
to  implement  high  performance  analog  computing  components. 

We  have  recently  begun  the  design  of  a  wide-dynamic  range  analog  multiplier 
based  on  pulse-frequency  encoding.  Specifically,  we  utilize  the  fact  that  the  amount 
of  charge  accumulated  by  a  CMOS  switched  capacitor  integrator  is  proportional 
to  the  product  of  an  analog  input  voltage  and  the  number  of  clock  cycles.  Thus 
by  utilizing  an  analog  voltage- to- frequency  converter  to  generate  a  clock  frequency 
that  is  directly  proportional  to  some  analog  input,  it  is  possible  to  implement  a 
two-quadrant  analog  multiplier  with  analog  inputs  and  analog  output.  The  ac¬ 
curacy  of  the  pulse-frequency  multiplier  will  be  limited  primarily  by  (i)  parasitic 
charge  injection  induced  by  the  switch  and  (ii)  the  maximum  integration  time  which 
imposes  a  minimum  clock  frequency.  A  full  four-quadrant  multiplier  can  be  imple¬ 
mented  by  using  a  more  sophisticated  CMOS  switched  capacitor  integrator  with 
additional  control  circuitry  to  represent  negative  values.  In  particular,  the  voltage- 
to-frequency  converter  needs  to  generate  a  clock  frequency  which  is  proportional  to 
the  absolute  magnitude  of  the  analog  input  and  also  provide  an  additional  signal 
indicating  the  sign  of  the  input  value.  Subsequently,  this  control  signal  can  be  used 
to  conditionally  reverse  polarity  of  the  input  charge. 

The  pulse-frequency  multiplier  is  well  suited  for  integration  into  a  massively 
parallel  processor  for  principal  component  analysis.  The  processing  components 
for  each  element  of  the  matrix  can  be  implemented  by  incorporating  multiple  in¬ 
put  structures  on  a  single  CMOS  switched  capacitor  integrator.  Thus  each  input 
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structure  provides  input  charge  at  a  rate  that  is  proportional  to  the  product  of  an 
analog  input  voltage  and  an  asynchronous  clock  frequency.  Furthermore,  the  matrix 
elements,  htJ  are  stored  as  the  charge  on  the  integration  capacitor  of  the  CMOS 
integrators.  By  proper  weighting  of  the  input  capacitors  in  each  input  structure,  it 
is  possible  to  implement  directly  the  computation  of  hi}  and  the  continued  update 
of  htj.  Component  block  and  transistor  level  simulations  of  the  pulse-frequency 
multiplier  and  the  massively  parallel  processor  axe  currently  being  performed. 
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m.6  Ordinal  Optimization.  Y.C.  Ho,  R.  Sreenivas,  Z.  Bo  Tang,  L.  Mol- 
lamustafaoglu,  Grants  N00014-89-J-1023,  N00014-89-J-0093,  NSF-CDR-88- 
03012,  NSF-ESC-91-02346,  DAAL-03-86-K-0171,  DAAL-03-91-G-0194;  Re¬ 
search  Unit  10. 

We  have  shown  that  the  order  of  performance  of  discrete  event  dynamic  sys¬ 
tems  (DEDS)  under  various  parameter  settings  are  relatively  immune  to  estimation 
errors.  Consequently,  if  we  concentrate  first  on  the  determination  of  a  subset  of 
good,  better,  and  best  designs  before  evaluating  the  actual  cardinal  performance 
values,  then  we  can  often  effect  a  great  increase  in  speed  of  the  optimization  search 
[1].  During  1991-92  we  further  developed  this  approach  with  respect  to 

(i)  Study  of  the  effect  of  correlation  on  ordinal  optimization  in  multiple  para¬ 
metrically  different  but  structurally  similar  simulations  [2]. 

(ii)  Study  of  the  use  of  adaptation  in  ordinal  optimization.  In  particular,  the  use 
of  genetic- like  algorithms  to  help  with  ordinal  search  was  investigated  [3]. 

Application  of  this  algorithm  has  been  made  to  the  optimal  design  of  robot 
arms.  The  PI  gave  a  lecture  course  in  Washington,  D.C.  to  personnel  from  the 
Army  Concepts  Analysis  Agency  and  the  Office  of  the  Joint  Chiefs  of  Staff  on 
5/29-30/92.  Possible  application  of  the  ordinal  optimization  algorithm  to  a  real 
logistics  problem  in  the  offices  of  JCS  and  to  the  wax  game  simulation  of  the  CAA 
are  under  discussion  as  a  result  of  the  lecture  course. 


III.7  Infinitesimal  Perturbation  Analysis  Y.C.  Ho,  Grant  ONR-N00014-89- 
J- 1023;  Research  Unit  10. 

A  unified  approach  to  IPA  was  completed.  We  now  have  a  very  complete 
story  of  the  computational  and  analytical  properties  of  IPA  in  the  sensitivity  cal¬ 
culation  of  DEDS  performances  [4-6]  integrating  the  recent  contributions  to  this 
topic  from  Russia  (Gaivironski)  and  France  (Bremaud). 
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IV.  ELECTROMAGNETIC  PHENOMENA 


Personnel 


Prof.  T.  T.  Wu 
Prof.  R.  W.  P.  King 
Dr.  J.  M.  Myers 
Dr.  H.-M.  Shen 


Mr.  G.  Fikioris 
Mr.  D.  K.  Freeman 
Mr.  V.  Houdzoumis 
Ms.  M.  Owens 
Ms.  B.  H.  Sandler 


Research  in  the  area  of  electromagnetic  radiation  is  directed  toward  the  so¬ 
lution  of  practical  problems  through  the  complete  understanding  of  the  underlying 
physical  phenomena.  This  involves  the  coordinated  application  of  modern  analyt¬ 
ical,  numerical,  and  experimental  techniques  and  the  use  of  high-speed  computers 
and  precision  instrumentation.  Application  is  also  made  of  modeling  techniques 
and  the  principle  of  similitude.  Most  practically  significant  problems  in  the  area 
are  sufficiently  complicated  that  extensive  computation  and  measurement  are  of¬ 
ten  required  to  justify  approximations  that  are  usually  necessary.  Where  possible, 
general  formulas  are  obtained  and  verified  experimentally  so  that  the  phenomenon 
under  study  can  be  understood  physically  in  analytical  form  and  not  just  as  a  set 
of  numbers. 

The  researches  are  concerned  primarily  with  the  properties  of  antennas  and 
arrays  and  of  the  electromagnetic  fields  they  generate  in  various  practically  im¬ 
portant  environments  that  lead  to  difficult  problems  with  complicated  boundary 
conditions.  Examples  include  the  properties  and  fields  of  dipoles,  insulated  anten¬ 
nas,  traveling-wave  antennas  and  arrays,  crossed  dipoles,  and  loops  when  located 
near  the  boundary  between  two  media  such  as  air  and  the  earth  or  sea,  or  the 
oceanic  crust  and  sea  water;  microstrip  patch  antennas  and  Beverage-type  anten¬ 
nas  for  over-the-horizon  radar;  and  resonant  circular  and  other  closed-loop  arrays  of 


67 


parallel  dipoles  in  free  space  and  on  the  earth  or  sea.  The  field  near  the  boundary 
includes  lateral  electromagnetic  waves  with  unusual  properties  that  require  detailed 
study  for  important  applications  in  remote  sensing  from  the  surface  of  the  earth 
or  sea,  the  arctic  ice  or  the  sea  floor,  and  in  horizontally-layered  media.  Also  of 
theoretical  and  practical  interest  are  the  generation,  propagation,  and  reception  of 
lateral  electromagnetic  pulses  along  boundaries;  and  solitary  electromagnetic  pulses 
with  slow  rates  of  decay. 

IV.  1  Electromagnetic  Field  of  Dipoles  and  Patch  Antennas  on  Microstrip. 

R.  W.  P.  King,  Grant  N00014-89-J-1023;  Research  Unit  11. 

Microstrip  antennas  consist  of  variously  shaped  and  excited  conducting  strips 
or  patches  on  the  dielectric  substrate.  Their  primary  function  is  to  generate  a 
significant  electromagnetic  field  at  distant  points.  This  field  is  the  superposition  of 
contributions  from  infinitesimal  elements  of  current  or  unit  electric  dipoles  in  the 
patch.  The  unit  dipole  has  an  electric  moment  Ih  =  1  Am;  when  multiplied  by 
I(x")dx",  the  dipole  is  an  element  of  a  patch  antenna.  The  complete  field  of  such 
an  antenna  is  obtained  by  an  integration  over  the  current  distribution.  This  has  not 
been  determined  analytically  for  commonly  used  shapes.  In  practice,  it  is  obtained 
by  numerical  methods  or  by  an  educated  guess.  Since  the  far  field  is  not  sensitive  to 
the  details  of  the  current  distribution,  a  usually  adequate  picture  of  the  radiating 
properties  of  a  patch  antenna  can  be  obtained  in  this  manner. 

Formulas  for  the  complete  electromagnetic  field  of  a  unit  horizontal  electric 
dipole  on  the  surface  of  a  dielectric-coated  conductor  (microstrip)  were  derived  last 
year  [1].  They  have  been  used  here  [2]  to  derive  expressions  for  the  complete  far 
field  of  an  x' -directed  element  at  an’’  point  p,<f>' ,z'  in  the  air;  this  includes  the  space 
wave  and  the  surface  wave.  From  these,  the  corresponding  formulas  for  the  far  field 
of  an  antiresonant  rectangular  patch  antenna  with  an  assumed  sinusoidal  current 
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distribution  have  been  derived. 


References: 

1.  R.  W.  P.  King,  “The  Electromagnetic  Field  of  a  Horizontal  Electric  Dipole  in 
the  Presence  of  a  Three-Layered  Region,”  J.  Appl.  Phys.  69(12),  7987-7995 
(1991). 

2.  R.  W.  P.  King,  “Electromagnetic  Field  of  Dipoles  and  Patch  Antennas  on 
Microstrip,”  Radio.  Sci.  27(1),  71-78  (1992). 

IV.2The  Circuit  Properties  and  Complete  Fields  of  Wave  Antennas  and 
Arrays.  R.  W.  P.  King,  Grant  N00014-89-J-1023;  Research  Unit  11. 

The  wave  antenna  is  important  for  communication  between  points  on  or  near 
the  surface  of  the  earth  and  especially  as  a  versatile  backscatter  radar  capable  of 
detecting  targets  in  the  whole  range  from  distant  and  high-flying  to  quite  close  and 
low-flying.  For  distant  targets  the  curvature  of  the  earth  makes  the  use  of  reflection 
from  the  ionosphere  a  necessary  complication.  Transmission  to  the  300  to  500  km 
distant  ionosphere  requires  quantitative  knowledge  of  the  field  at  all  upward  angles. 
When  the  target  is  not  an  ICBM  thousands  of  kilometers  away  but  a  cruise-missile 
launched  by  a  submarine  or  a  low-flying  drug-carrying  aircraft  within  1000  km  of  the 
radar  transmitter-receiver,  reflections  from  the  ionosphere  are  best  supplemented 
by  the  surface  wave.  This  is  unaffected  by  geomagnetic  storms  and  independent  of 
skip  distances. 

The  wave  antenna  consists  of  a  long  horizontal  wire  in  air  electrically  very 
close  to  the  surface  of  the  earth  or  sea  and  so  terminated  that  a  pure  traveling 
wave  propagates  along  it  in  both  the  transmitting  and  receiving  modes.  This  can 
be  accomplished  with  suitable  ground  connections  in  series  with  a  resistor  as  in  the 
familiar  Beverage  antenna.  Alternatively,  each  termination  may  be  a  resistor  in 
series  with  a  horizontal,  quarter-wave  monopole  to  form  a  horizontal-wire  antenna. 
With  either  form,  the  generator-receiver  is  located  adjacent  to  the  resistor  at  one 
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end.  The  wave  antenna  was  analyzed  previously  [1]  subject  to  the  condition  p2 
z'2  where  p  is  the  radial  distance  from  the  source  and  z'  =  —z  is  the  height  of  the 
point  of  observation  above  the  surface  of  the  earth. 

In  a  recent  investigation  [2],  the  complete  electromagnetic  field  generated  by 
a  wave  antenna — including  both  the  space  wave  and  the  surface  wave  at  all  points 
in  the  air  and  on  the  surface  of  the  earth  or  sea — has  been  derived.  A  knowledge  of 
this  complete  field  is  essential  for  an  over-the-horizon  radar  that  makes  use  of  both 
ionospheric  reflection  for  detecting  distant  high-flying  targets  and  the  surface  wave 
for  detecting  near  and  low-flying  targets. 

Although  the  electromagnetic  field  of  a  unit  electric  dipole  is  smaller  by  the 
factor  |  k2/kl  |  when  it  is  horizontal  at  a  small  height  d  than  when  it  is  vertical,  a  long 
traveling-wave  antenna  with  a  very  substantial  directive  gain  that  can  approximate 
or  exceed  the  ratio  \k1/k2\  can  be  constructed  of  horizontal  elements.  Indeed,  an 
array  of  broadside  and  collinear  wave  antennas  can  provide  an  electromagnetic 
field  of  much  greater  intensity  and  at  much  lower  constructional  cost  than  can 
be  achieved  with  vertical  monopoles  alone.  This  is  not  obvious  from  the  earlier 
analysis  [1]  since  it  determined  only  the  surface-wave  field  of  a  single  Beverage  or 
horizontal-wire  antenna. 

It  is  also  possible  to  combine  a  radial  array  of  individually  directive  wave 
antennas  with  a  vertical  monopole  to  obtain  a  useful  omnidirectional  antenna  that 
generates  a  powerful  surface  wave  combined  with  a  significant  space  wave.  In  [2], 
the  design  of  such  an  array  is  presented.  It  consists  of  a  vertical  monopole  and 
IV  =  10  horizontal  wave  antennas  at  a  small  height  d  =  45  cm  over  the  earth.  No 
ground  connections  axe  required.  It  is  shown  that  the  use  of  am  array  of  hor’  ontal 
wave  antennas  with  a  vertical  monopole  in  place  of  the  usual  grounding  network 
not  only  increases  the  space  wave  but  uses  the  power  ultimately  dissipated  in  the 
earth  to  generate  a  large  and  useful  surface  wave. 
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In  some  applications  in  radar  and  in  point-to-point  communication,  a  directive 
array  is- preferred  over  an  omnidirectional  one.  In  such  cases,  an  antenna  that  radi¬ 
ates  a  large  field  within  a  specified  angle  and  a  small  field  in  ail  other  directions  is 
advantageous.  Such  an  array  is  readily  obtained  by  a  simple  modification  of  the  om¬ 
nidirectional  array  discussed  above.  Other  combinations  of  vertical  monopoles  and 
horizontal  wave  antennas  have  interesting  and  useful  characteristics.  An  example  is 
a  highly  directive  broadside  array  in  which  each  element  is  a  single  vertical  antires¬ 
onant  monopole  in  series  with  a  single  horizontal  wave  antenna.  A  base-insulated 
resonant  dipole  is  used  as  a  parasitic  reflector  for  each  vertical  element. 
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IV.3The  Propagation  of  a  Gaussian  Pulse  in  Sea  Water.  R.  W.  P.  King, 
B.  H.  Sandler,  and  M.  Owens,  Grant  N00014-89-J-1023  and  Raytheon  Co.  (un¬ 
der  Contract  N61533-90-C-0080  with  ONR,  Annapolis  Detachment);  Research 
Unit  11. 

The  propagation  of  electromagnetic  waves  in  sea  water  is  much  more  compli¬ 
cated  than  in  lake  water  because  the  attenuation  and  phase  constants  are  practically 
equal  and  both  depend  strongly  on  the  frequency  in  a  nonlinear  manner.  The  in¬ 
teresting  question  arises:  How  does  a  single  pulse  progress  in  sea  water?  Is  such 
a  pulse  a  useful  vehicle  for  remote  sensing  in  the  ocean?  These  questions  can  be 
investigated  conveniently  with  a  Gaussian  pulse  of  current  applied  to  a  horizon¬ 
tal  electric  dipole  on  the  surface  of  the  sea.  The  Gaussian  pulse  is  essentially  a 
low-frequency  pulse  in  the  sense  that  the  amplitudes  of  the  constituent  frequencies 
decrease  very  rapidly  with  increasing  frequency.  It  is  clearly  well-suited  to  propaga- 
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tion  in  sea  water  where  the  exponential  attenuation  decreases  with  the  frequency. 
When  a  Gaussian  pulse  is  applied  to  a  horizontal  electric  dipole  on  the  surface  of 
the  sea,  aui  electromagnetic  pulse  is  generated  that  propagates  upward  into  the  air, 
downward  into  the  sea,  and  horizontally  along  the  boundary  in  a  surface  wave.  Of 
interest  here  is  the  pulse  propagating  vertically  downward  and  impinging,  for  ex¬ 
ample,  on  a  horizontal  metal  cylinder.  Can  the  cylinder  be  detected  by  observing 
the  backscattered  pulse? 

Detailed  quantitative  information  has  been  obtained  for  the  properties  of  a 
Gaussian  pulse  propagating  in  sea  water  [1].  Of  particular  significance  are  the  fol¬ 
lowing  points:  1)  The  amplitude  of  the  propagating  pulse  decreases  very  rapidly 
with  distance.  2)  The  shape  of  the  pulse  changes  continuously  with  distance  from 
the  initial  Gaussian  shape  to  that  of  the  time-derivative  of  the  Gaussian.  3)  The  ve¬ 
locity  of  propagation  of  electromagnetic  waves  in  sea  water  is  frequency-dependent, 
decreasing  with  frequency.  As  the  pulse  propagates,  the  higher  frequencies  in  its 
spectrum  are  more  rapidly  attenuated  than  the  lower  frequencies  so  that  the  higher- 
frequency  content  of  the  pulse  decreases  rapidly  and  with  it  the  velocity  of  prop¬ 
agation  of  the  pulse.  4)  The  apparent  velocity  of  propagation  of  the  maximum 
magnitude  of  the  Gaussian  pulse  is  anomalous.  Because  of  the  change  in  shape 
from  the  Gaussian  to  that  of  the  derivative,  the  maximum  of  the  pulse  moves  back¬ 
wards  for  a  short  time.  Actually  all  components  move  forward,  only  the  maximum 
of  their  sum  moves  backwards. 

To  study  the  possible  use  of  such  a  pulse  in  the  detection  of  scattering  ob¬ 
stacles,  a  direct  transmission-reflection  method  without  the  intermediary  of  surface 
waves  has  been  used.  This  makes  use  of  a  horizontal  electric  dipole  in  the  form  of 
an  insulated  conductor  with  the  length  2 h  terminated  in  bare  ends.  The  dipole  is 
electrically  short  for  all  frequencies  involved  in  the  Gaussian  current  pulse  that  is 
generated  in  it.  The  electromagnetic-field  pulse  travels  vertically  down  into  the  sea. 
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It  induces  a  current  in  a  metal  cylinder  with  the  length  100  m  and  cross-sectional 
diameter  12.5  m  located  at  a  depth  z  —  100  m.  The  current  in  the  cylinder  radiates 
an  electromagnetic  field  in  all  directions  in  the  sea  including  upward.  The  signal 
reaching  a  receiver  adjacent  to  the  transmitting  dipole  is  the  means  for  detecting 
the  cylinder. 

The  shape  and  amplitude  of  the  pulse  reflected  by  the  cylinder  depends  on  the 
spectrum  of  frequencies  involved  and,  therefore,  on  the  width  of  the  incident  pulse. 
The  current  induced  in  the  cylinder  by  each  component  of  frequency  is  given  in  [2]. 
In  the  general  case,  the  field  scattered  by  each  component  of  frequency  has  to  be 
determined  as  in  [1]  and  then  summed  over  all  frequencies.  In  a  preliminary  study 
this  step  can  be  avoided  if  the  range  of  relevant  frequencies  is  restricted  to  those  for 
which  the  length  of  the  cylinder  is  electrically  short  in  sea  water.  Gaussian  pulses 
with  widths  2 tx  >0.1  sec  satisfy  this  condition.  For  them,  the  current  induced 
in  the  cylinder  has  the  same  frequency  dependence  as  the  incident  field  so  that 
the  scattered  field  will  also  have  this  dependence.  This  means  that  an  incident 
Gaussian-field  pulse  should  be  scattered  as  a  Gaussian- field  pulse. 

If  significant  higher  frequencies  are  present,  the  induced  current  becomes 
frequency- dependent  and  the  scattered  field  is  modified.  Actually  the  induced  cur¬ 
rent  distribution  differs  little  from  the  parabolic  shape  for  electrically  much  longer 
scattering  elements.  The  principal  difference  is  an  increase  in  the  amplitude  of  the 
induced  current  as  the  scatterer  approaches  the  resonant  length  near  0xh  =  x/2. 
It  follows  that  the  range  of  pulse  widths  can  be  decreased  to  2tx  =  0.01  sec  with 
the  understanding  that  a  scattered  pulse  will  preserve  the  shape  of  the  incident- 
field  pulse  quite  accurately,  but  that  the  actual  scattered-field  amplitude  should  be 
somewhat  greater  than  that  calculated  using  the  approximate  current. 

It  has  been  found  that  a  cylinder  at  a  depth  of  100  m  is  detectable  with  a 
Gaussian  pulse  of  width  2tt  =  0.01  sec,  but  not  without  a  complicated  cancellation 
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of  the  large  direct  field  when  a  Gaussian  pulse  of  width  2tl  =  1  sec  is  used.  A 
complete  investigation  of  a  range  of  relevant  depths  and  pulse  widths  can  be  made. 
However,  for  pulses  shorter  than  2tj  =  0.01  sec,  the  general  formula  for  the  current 
induced  in  the  cylinder  must  be  used  and,  since  this  is  frequency- dependent,  the 
response  to  an  incident  pulse  requires  numerical  evaluation. 
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IV.4The  Propagation  of  a  Low-Frequency  Burst  in  Sea  Water.  R.  W.  P. 
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Research  Unit  11. 

The  study  of  the  propagation  of  single  Gaussian  pulses  in  sea  water  [1] — 
discussed  under  Topic  IV. 3 — has  been  extended  to  include  the  propagation  in  sea 
water  of  a  short  burst  of  oscillations  when  its  amplitude  is  limited  by  a  Gaussian 
envelope.  In  particular,  a  complete  analysis  has  been  made  [2]  of  the  propagation 
in  sea  water  of  a  burst  of  25  cycles  at  /  ~  25  Hz  with  amplitudes  limited  by  a 
Gaussian  envelope  with  the  half-width  tx  =  0.5  sec.  The  properties  of  the  burst  in 
transmission  are  governed  in  part  by  the  Gaussian  envelope,  in  part  by  the  burst 
of  oscillations.  The  most  important  characteristic  is  the  very  rapid  decrease  in  the 
overall  amplitude  due  to  the  electrically  short  distances  involved  and  the  initial  1/z3 
dependence  on  distance.  Owing  to  the  very  low  frequencies  involved  in  the  burst 
and  the  spectrum  of  the  envelope,  the  exponential  attenuation  is  not  great. 

Just  as  for  the  isolated  Gaussian  pulse  [1],  the  apparent  continuity  of  Gaussian 
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shape  is  misleading.  Actually,  the  initial  burst  begins  with  a  Gaussian  envelope  but 
the  entire  configuration  changes  rapidly  to  its  time-derivative.  In  the  process  the 
maximum  of  the  burst  shifts  to  the  initial  maximum  of  the  time-derivative.  The 
backward  movement  of  the  maximum  can  be  observed  but  this  is  much  more  diffi¬ 
cult  owing  to  the  oscillations;  also  it  is  much  more  restricted  to  very  short  distances. 
A  careful  study  of  the  numerical  evaluation  of  the  amplitudes  of  the  individual  os¬ 
cillations  shows  that  the  apparent  velocity  of  propagation  va  =  z/t  of  the  maximum 
of  the  envelope  decreases  from  large  values  near  z'  =  z  =  0.1  to  an 

approximately  constant  value  of  the  order  of  magnitude  of  the  phase  velocity  for  a 
steady-state  traveling  wave  in  sea  water  at  /  ~  25  Hz.  This  behavior  differs  from 
that  of  the  single  Gaussian  pulse  which  continuously  loses  higher  frequencies  from 
its  spectrum  as  it  travels  in  sea  water  and,  hence,  advances  more  and  more  slowly. 
The  burst  has  primarily  the  properties  of  the  25  Hz  oscillations  at  all  distances. 

Primary  interest  is  in  the  electric  field  observed  at  a  receiver  on  the  air-sea 
surface  near  the  transmitting  dipole.  As  for  the  single  Gaussian  pulse,  it  is  assumed 
that  the  back-scattered  field  is  approximately  the  field  incident  on  the  scatterer 
but  reversed  in  direction.  The  scattered  burst  is  smaller  than  the  direct  burst 
observed  at  the  receiver  on  the  surface  by  a  factor  of  the  order  10~4.  It  arrives  at 
the  receiver  while  the  direct  burst  is  still  rising  to  its  maximum.  If  the  direct  burst 
can  be  balanced  out  in  its  entirety  at  the  receiver,  the  scattered  burst  is  sufficiently 
large  to  be  detectable. 
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IV.5  The  Propagation  of  a  Radar  Pulse  in  Sea  Water.  R.  W.  P.  King,  T.  T. 
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Research  Unit  11. 

Conventional  radar  transmits  and  receives  wave  packets  that  approximate 
finite  sections  of  a  continuous-wave  field.  In  air,  these  propagate  with  no  change 
in  shape  since  all  frequencies  in  the  spectrum  of  the  modulated  pulse  travel  with 
the  velocity  of  light.  When  a  similar  wave  packet  is  transmitted  in  sea  water,  a 
much  more  complicated  phenomenon  occurs  because  the  medium  is  dissipative  and 
dispersive  so  that  the  wave  number  is  complex  and  not  linear  in  the  frequency.  As 
a  consequence,  the  rise  time  and,  with  it,  the  shape  of  the  wave  packet  are  greatly 
modified  as  the  pulse  travels  in  the  attenuating  and  dispersive  salt  water.  In  a 
new  study,  the  propagation  of  wave  packets  has  been  investigated  when  the  carrier 
frequency  is  sufficiently  low  to  permit  a  significant  distance  of  travel  before  the 
amplitude  becomes  too  small  to  provide  a  detectable  signal. 

As  a  necessary  preliminary  to  the  determination  of  the  electromagnetic  field 
scattered  by  a  metal  cylinder  submerged  in  the  ocean,  the  propagation  of  a  low- 
frequency  pulse  with  a  rectangular  envelope  has  been  evaluated  numerically  as  a 
function  of  the  distance  of  travel  [lj.  The  electric  fields  generated  by  an  electrically 
short  dipole  in  sea  water  excited  by  a  wave  packet  of  25.5  cycles  at  25.5  Hz  reveal 
a  complicated  behavior  as  the  packet  penetrates  to  increasing  depths.  The  shape 
of  the  wave  packet  is  increasingly  altered  with  distance  of  travel.  It  follows  that 
the  received  signal  is  quite  different  from  the  transmitted  signal  when  distances  are 
large. 

To  understand  the  numerically  evaluated  results,  an  analytical  investigation 
to  explain  the  computed  properties  was  carried  out.  Specifically,  the  propagation 
in  sea  water  of  an  electromagnetic  field  in  the  form  of  a  semi-infinite  wave  train 
or  a  radar  pulse  generated  by  an  electric  dipole  has  been  investigated  analytically 
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for  low  frequencies.  In  the  resulting  paper  [2],  the  frequency-domain  formula  for 
the  downward-traveling  field  of  a  horizontal  electric  dipole  excited  by  a  sinusoidally 
modulated  electric-current  pulse  is  Fourier  transformed  to  obtain  an  explicit  expres¬ 
sion  for  the  field  at  any  distance  in  the  time  domain.  Specific  application  is  made 
to  a  wave  packet  of  25-5  cycles  in  a  time  duration  of  one  second.  The  amplitude 
and  phase  velocity  of  the  wave  packet  are  determined  together  with  the  amplitudes 
of  the  initial  and  final  precursors.  Graphs  axe  displayed  and  discussed  for  a  range 
of  distances;  these  show  that  the  amplitude  of  the  wave  packet  decays  more  rapidly 
than  the  amplitudes  of  the  precursors. 

Supplementing  the  recent  study  by  Oughston  [3]  which  obtains  “an  asymptotic 
description  for  large  values  of  the  propagation  distance,”  the  present  study  provides 
a  complete  formula  for  the  outward-traveling  field  both  near  and  far  in  the  equatorial 
plane  of  a  dipole  source.  For  possible  application  to  the  detection  of  targets  in  the 
ocean  with  conventional  radar  techniques  instead  of  lateral  waves  as  described  in 
[4],  the  carrier  frequency  /0,  the  pulse  length  2tr,  and  the  repetition  rate  must  be 
selected  with  due  consideration  of  the  signal  velocity  and  attenuation  of  the  wave 
packet  and  of  its  distortion  due  to  the  leading-  and  trailing-edge  precursors.  These 
latter  are  attenuated  with  distance  much  more  slowly  than  the  wave  packet  and 
ultimately  dominate  the  advancing  pulse. 
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IV.6  Radiation  of  Electromagnetic  Waves  from  a  Tether  in  the  Iono¬ 
sphere  to  the  Surface  of  the  Earth.  R.  W.  P.  King,  Grant  N00014- 
89-J-1023  and  Raytheon  Co.  (under  Contract  N61533-90-C-0080  with  ONR, 
Annapolis  Detachment);  Research  Unit  11. 

The  possibility  of  transmitting  electromagnetic  waves  from  an  antenna  at¬ 
tached  to  an  orbiting  space  shuttle  to  the  surface  of  the  earth  has  been  suggested. 
The  antenna  generating  the  electromagnetic  field  would  consist  of  a  vertical  wire 
(the  tether),  driven  with  the  space  shuttle  as  ground.  The  electromagnetic  aspects 
of  the  problem  have  been  examined  with  available  knowledge  as  the  basis.  Because 
the  existing  theory  of  the  properties  of  the  ionosphere  and  of  antennas  moving  in  it 
is  approximate  and  complicated,  quantitative  results  can  be  obtained  only  in  terms 
of  a  relatively  simple  model.  The  specific  problem  investigated  is  the  electromag¬ 
netic  field  on  the  surface  of  the  earth  generated  by  currents  in  a  vertical  antenna 
moving  in  the  F-layer  of  the  ionosphere  at  a  height  of  400  km.  The  length  of  the 
antenna  is  4  km.  It  is  driven  at  its  upper  end  by  a  generator  voltage  VQ  against  the 
space  shuttle  as  a  ground.  The  ionosphere  is  assumed  to  extend  from  a  height  of 
150  km  to  infinity  as  a  homogeneous  medium.  It  is  given  a  sharp  boundary  with 
air  as  a  simplification  of  the  gradual  layered  boundary  between  100  and  200  km. 

The  properties  of  the  ionosphere  and  three  types  of  propagating  waves  have 
been  examined  for  two  frequency  ranges:  VLF  (9  <  /  <  30  kHz)  and  ELF  (30  < 
/  <  300  Hz).  Electromagnetic  plasma  waves  and  the  Whistler  mode  propagate  in 
the  entire  VLF  band,  and  Alfven  waves  in  the  ELF  band.  With  the  help  of  available 
approximate  formulas,  it  is  shown  that  a  measurable  incident  electric  field  can  be 
maintained  on  the  surface  of  the  earth. 
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IV.7  Electromagnetic  Fields  of  Vertical  and  Horizontal  Magnetic  Dipoles 
in  and  over  the  Sea.  R.  W.  P.  King  and  M.  Owens,  Grant  N00014-89-J- 
1023  and  Raytheon  Co.  (under  Contract  N61533-90-C-0080  with  ONR,  An¬ 
napolis  Detachment);  Research  Unit  11. 

Owing  to  the  symmetry  of  Maxwell’s  equations  in  free  space,  the  complete 
electromagnetic  field  of  the  magnetic  dipole  can  be  obtained  from  that  of  the  elec¬ 
tric  dipole  with  the  simple  substitutions:  E'  =  B//i0  and  B'  =  — eE,  where  E 
and  B  are  the  fields  of  the  electric  dipole,  E'  and  B'  the  fields  of  the  magnetic 
dipole.  When  an  electric  or  magnetic  dipole  is  near  the  plane  boundary  between 
two  electrically  different  regions,  the  electromagnetic  fields  generated  by  the  electric 
and  magnetic  dipoles  are  no  longer  simply  related  by  these  complementary  formu¬ 
las.  This  is  a  consequence  of  the  fact  that  the  boundary  conditions  for  the  electric 
and  magnetic  fields  on  the  boundary  between  the  two  different  regions  are  not  the 
same.  All  tangential  components  are  continuous,  but  the  normal  component  of  the 
electric  field  is  discontinuous  while  the  normal  component  of  the  magnetic  field  is 
continuous. 

A  complete  and  detailed  analysis  has  been  made  of  the  electromagnetic  fields 
generated  by  horizontal  and  vertical  electric  dipoles  near  a  plane  boundary  between 
electrically  different  regions  such  as  air  and  sea  water  [1]  (See  Topic  IV. 9).  The 
derivations  are  complex  and  elaborate.  In  order  to  obtain  the  corresponding  for¬ 
mulas  for  the  fields  of  magnetic  dipoles,  two  alternatives  are  available.  One  is  to 
parallel  the  procedures  carried  out  for  the  electric  dipoles.  The  other  is  to  combine 
the  fields  of  electric  dipoles  to  obtain  the  field  of  a  small  square  constructed  of  elec¬ 
tric  dipoles,  each  with  the  electric  moment  17,  and  to  replace  this  by  the  equivalent 
magnetic  moment  VI  =  (iu>n0l)  x  (II)  =  iu^0l2I.  Neither  procedure  is  simple,  but 
the  latter  is  less  complicated  in  that  it  involves  no  new  integrals  that  need  to  be 
evaluated. 

This  latter  approach  is  presently  being  carried  out.  The  y-directed  magnetic 
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dipole  is  represented  by  a  square  loop  of  electric  dipoles;  its  field  is  the  vector  sum 
of  the  fields  of  the  four  dipoles.  In  order  to  explore  the  procedure  and  check  its 
accuracy,  it  has  been  applied  first  to  the  fields  of  isolated  electric  dipoles,  and  the 
corresponding  known  field  of  the  magnetic  dipole  has  been  successfully  derived. 
Work  is  currently  underway  to  derive  the  six  components  of  the  field  of  a  y-directed 
horizontal  magnetic  dipole  with  its  center  at  a  height  d  in  air  over  the  surface  of 
the  earth  or  sea. 
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IV.8  The  Detection  of  Dielectric  Spheres  Submerged  in  Water.  R.  W.  P. 
King  and  S.  S.  Sandler,  Grant  N00014-89-J-1023  and  the  Small  Business  Ad¬ 
ministration  under  the  SBIR  Program  in  collaboration  with  Geo-Centers,  Inc., 
and  the  U.S.  Army  Medical  Branch;  Research  Unit  11. 

The  use  of  lateral  electromagnetic  waves  to  locate  buried  or  submerged  metal 
objects  in  the  earth  or  sea  with  the  horizontal  electric  dipole  as  the  source  has 
been  described  [l]-[3].  A  related  problem  is  the  corresponding  search  for  dielectric 
objects  embedded  in  a  general  dissipative  medium.  Examples  are  plastic  mines  in 
the  ocean  and  fragments  of  exploded  plastic  land  mines  embedded  in  human  flesh. 
The  important  difference  between  the  remote  sensing  for  metal  and  plastic  targets 
is  the  inverse  nature  of  the  relative  magnitudes  of  the  wave  numbers  characteristic 
of  electromagnetic  wave  propagation  in  the  ambient  medium  and  in  the  embedded 
target.  For  simplicity  the  target  can  be  taken  to  be  a  sphere  consisting  of  air  bubbles 
or  plastic  embedded  in  lake  water,  sea  water  or  human  flesh.  The  purpose  of  the 
study  is  to  investigate  the  possibility  of  locating  dielectric  spheres  in  water  or  flesh 
by  means  of  observations  made  on  the  boundary  surface  in  the  air. 
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The  source  is  a  horizontal  electric  dipole  in  the  air  on  its  boundary  with  the 
ambient  medium  that  contains  the  target.  It  generates  an  electromagnetic  field 
that  propagates  in  the  air  along  the  boundary  as  a  surface  wave  of  a  type  known 
as  a  lateral  wave.  The  amplitude  fronts  of  the  field  in  the  earth  or  water  propagate 
vertically  downward;  the  phase  fronts  travel  at  an  angle.  It  is  this  field  that  is 
incident  on  a  target  located  at  a  depth  z.  In  practice  the  field  is  not  maintained 
by  an  infinitesimal  unit  dipole,  but  by  a  center-driven  antenna  located  parallel  to 
and  at  a  small  height  in  the  air  above  the  plane  boundary.  The  source  is,  in  effect, 
a  center-driven  eccentrically  insulated  dipole.  If  the  dipole  antenna  is  adjusted  in 
length  to  be  resonant  and  the  height  is  carefully  chosen,  the  driving-point  impedance 
will  be  resistive  and,  hence,  readily  matched  to  a  feeding  transmission  line. 

The  basic  theory  has  been  developed  [4]  for  the  use  of  horizontal  electric 
dipoles  on  the  surface  of  a  dissipative  region  to  detect  small  dielectric  spheres  em¬ 
bedded  at  some  depth  in  that  region.  Since  interest  is  in  the  detection  of  quite 
small  targets  and  the  scattering  by  dielectrics  is  much  weaker  than  by  metals,  it  is 
necessary  to  make  use  of  quite  high  frequencies  in  order  to  have  adequate  resolving 
power.  It  has  been  shown  that  detectable  reflected  signals  can  be  generated  with 
acceptable  power  so  that  it  should  be  possible  to  locate  small  dielectric  objects  with 
suitably  designed  apparatus. 
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4.  R.  W.  P.  King  and  S.  S.  Sandier,  “The  Detection  of  Dielectric  Spheres  Sub¬ 
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IV.9  Lateral  Electromagnetic  Waves  (book).  R.  W.  P.  King,  M.  Owens,  and 
T.  T.  Wu,  Grant  N00014-89-J-1023;  Research  Unit  11. 

A  new  746-page  book  entitled  Lateral  Electromagnetic  Waves:  Theory  and 
Applications  to  Communications,  Geophysical  Exploration,  and  Remote  Sensing, 
by  R.  W.  P.  King,  M.  Owens,  and  T.  T.  Wu  has  been  published  [1].  This  book 
analyzes  electromagnetic  radiation  due  to  sources  located  near  a  boundary  between 
two  media  like  air  and  earth  or  air  and  sea.  Such  radiation  has  properties  that  axe 
different  from  those  familiar  from  the  study  of  incident  plane  waves.  Of  these,  the 
most  important  property  is  a  surface  wave  of  a  type  known  as  a  lateral  wave,  which 
propagates  along  the  boundary.  (Analogous  waves  that  occur  in  seismology  are 
called  head  waves.)  Electromagnetic  lateral  waves  have  many  applications  which 
are  treated  in  this  book.  They  include  the  investigation  of  the  properties  of  the 
Earth’s  mantle,  radio  communication  with  submarines  and  over  the  Earth’s  surface, 
remote  sensing  and  over-the-horizon  radar.  More  general  situations,  in  which  there 
axe  more  than  two  layers,  axe  also  treated.  These  include  the  effect  of  sediment  in 
the  study  of  the  oceanic  crust  and  of  the  Arctic  ice  in  the  detection  of  submerged 
submarines.  The  analysis  of  microstrip  circuits  and  antennas  involves  three  layers, 
namely,  the  air,  the  dielectric  substrate,  and  the  conducting  base. 

The  systematic  study  presented  in  this  book  begins  with  the  rigorous  analyt¬ 
ical  derivation  of  the  complete  set  of  electromagnetic  waves — spherical  waves  and 
lateral  waves — generated  by  horizontal  and  vertical  dipoles  near  or  on  the  plane 
boundary  between  two  half-spaces  with  quite  different  characteristic  wave  veloci¬ 
ties.  Formulas  are  first  derived  in  the  form  of  general  integrals,  which  constitute 
the  exact  solution.  From  these,  simple  continuous  formulas  are  obtained  that  are 
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accurate  and  valid  everywhere  except  very  close  to  the  source.  These  simple  and 
practical  formulas  are  available  in  a  book  for  the  first  time. 

The  book  contains  the  results  of  current  and  past  researches  on  lateral  elec¬ 
tromagnetic  waves  that  have  been  supported  in  large  part  by  the  Joint  Services 
Electronics  Program.  It  is  characterized  by  its  broad  generality  and  mathematical 
exactitude,  and  by  its  detailed  treatment  of  important  practical  applications.  It 
should  be  a  very  useful  tool  and  comprehensive  reference  for  researchers  in  atmo¬ 
spheric  physics,  geophysics,  electrical  engineering,  applied  mathematics,  and  math¬ 
ematical  physics. 
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IV.10  Theoretical  Study  of  Electromagnetic  Pulses  with  a  Slow  Rate  of 
Decay.  T.  T.  Wu,  J.  M.  Myers,  H.  M.  Shen,  R.  W.  P.  King,  and  V.  Houd- 
zoumis,  Grant  N00014-89-J-1023  and  Army  LABCOM  Contract  DAAL02- 
89-K-0097;  Research  Unit  11. 

The  interesting  possibility  of  generating  electromagnetic  pulses  of  finite  total 
radiated  energy  that  decreases  with  distance  much  more  slowly  than  the  usual  r~ 2 
has  been  investigated  both  theoretically  (this  topic)  and  experimentally  (see  Topic 
IV.  11).  Such  pulses  are  conveniently  referred  to  as  electromagnetic  (EM)  missiles 
[1].  Theoretical  work  during  the  past  year  has  been  concerned  with  the  scattering 
of  EM  missiles  from  various  objects.  A  paper  has  been  published  [2]  that  analyzes 
the  two-dimensional  problem  of  the  scattered  energy  flux  due  to  an  EM  missile 
incident  upon  an  infinite  circular  perfectly  conducting  cylinder.  The  analysis  is 
based  on  an  expansion  of  the  field  in  terms  of  Bessel  functions.  The  paper  also 
treats  the  incident  missile  locally  as  a  plane-wave  transient  impinging  on  the  point 
of  specular  reflection,  and  applies  geometrical  optics  to  compute  the  reflection.  This 
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“geometric  approximation”  is  found  to  agree  exactly  with  the  explicit  calculation 
of  fields,  thus  justifying  both  the  local  approximation  by  a  plane  wave  and  the 
use  of  geometrical  optics.  The  geometrical  method  is  generalizable  to  a  broad 
class  of  convex  objects  with  nonzero,  bounded  surface  curvature.  The  result  for 
the  sphere  is  that  the  backscattering  preserves  linear  polarization,  reverses  circular 
polarization,  and  returns  an  energy  flux  proportional  to  r~(2+<),  where  the  incident 
missile  delivers  a  flux  r~e.  The  backscattered  energy  varies  as  the  square  of  the 
radius  of  the  sphere. 

The  r-2-<  behavior  of  the  energy  of  EM  missiles  scattered  from  spheres  differs 
drastically  from  the  r-t  dependence  of  the  backscattered  energy  on  range  from  a 
flat  plate.  The  question  arises  whether  there  are  objects  that  exhibit  a  scattering 
behavior  that  is  intermediate  between  that  for  round  and  flat  objects.  A  theoretical 
study  has  shown  that  such  objects  do  exist.  To  study  their  properties,  it  is  conve¬ 
nient  to  start  in  two  dimensions,  where  round  objects  echo  as  r”1"*,  in  place  of  the 
r~2~e  for  three  dimensions.  In  two  dimensions,  the  generalized  super-ellipse  defined 
as  |x|l/  +  |y|"  =  1  is  such  an  object  when  v  >  2.  The  range  dependence  of  the  energy 
backscattered  by  a  super-ellipse  illuminated  by  an  EM  missile  has  been  derived  and 
shown  to  be  intermediate  between  the  two  previously  reported  cases.  An  asymptotic 
formula  for  the  backscattered  energy  has  also  been  derived.  This  formula  is  valid 
for  large  distances  between  the  source  and  scatterer  and  is  uniform  in  e.  It  shows 
that  the  qualitative  behavior  of  the  backscattered  energy  depends  on  the  quantity 
q  =  1  —  2v~l  -  e.  In  particular,  for  q  >  0,  the  energy  attenuates  as 
whereas  for  q  <  0  it  attenuates  as  r_1+*.  When  q  =  0,  the  energy  attenuates  as 
r-1  lnr.  The  behavior  in  the  transitional  region  near  q  —  0  can  also  be  obtained 
from  the  asymptotic  formula.  Similar  results  can  be  shown  for  three-dimensional 
super-ellipsoids,  where  the  characteristic  quantity  is  then  q  —  2  —  4i/-1  —  e. 
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IV.ll  Experimental  Study  of  Electromagnetic  Pulses  with  a  Slow  Rate 
of  Decay.  H.-M.  Shen,  R.  W.  P.  King,  and  T.  T.  Wu,  Grant  N00014-89- 
J-1023  and  Army  LABCOM  Contract  DAAL02-89-K-0097;  Research  Unit 
11. 

An  improved  L-antenna  has  been  designed  to  measure  a  long  pulse  or  CW  elec¬ 
tromagnetic  field.  The  new  design  makes  use  of  a  distributed  resistance  to  eliminate 
the  undesirable  multiple  parasitic  reflections.  Two  methods  have  been  employed — 
one  more  accurate  than  the  other — to  analyze  a  two-wire  parallel  transmission  line 
with  distributed  resistance  and,  based  on  the  results  of  these  theoretical  models,  L- 
antennas  with  a  distributed  resistance  have  been  constructed.  Measurements  made 
with  the  improved  L- antennas  clearly  show  that  the  presence  of  the  distributed 
resistance  significantly  reduces  the  reflected  pulses.  The  design  based  on  the  more 
accurate  of  the  two  methods  produces  the  better  results. 

Also  during  this  reporting  period,  laser-triggered  electrical  discharges  have 
been  reviewed  as  a  mechanism  for  generating  EM  missiles  [1],  Since  1988,  Auston, 
Zhang,  and  their  coworkers  at  Columbia  University  [2]  have  used  laser-triggering  to 
provide  electromagnetic  current  over  an  element  of  area  on  the  order  of  1  cm2,  with 
demonstrated  subpicosecond  electromagnetic  pulses  free  of  any  carrier,  some  with 
bandwidths  of  1012  Hz,  and  some,  of  less  bandwidth,  with  peak  power  densities  of 
64  W/cm2.  These  advances  have  been  outlined  in  [1]  and  some  of  the  problems 
and  prospects  for  coordinating  many  such  elements  are  discussed.  A  design  for 
a  launcher  of  electromagnetic  missiles  for  a  target  100  km  distant  is  discussed  in 
terms  of  the  time  dependence  of  a  current  pulse  over  an  aperture.  Two  issues  in 
implementing  this  design  are:  1 )  to  provide  the  required  current  distribution  over  a 
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small  surface  element;  and  2)  to  coordinate  many  such  elements  in  an  array.  A  novel 
approach  for  coordinating  pulses,  based  on  an  adaptive  technique  already  used  for 
pulsed  lasers,  is  proposed  as  a  means  for  measuring  and  correcting  critical  timing 
discrepancies  in  an  array  of  pulse  generators;  beam  focusing  can  also  be  controlled 
in  this  way. 
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IV.12  Lessons  for  Pulsed-Array  Radar  from  Quantum  Light  Detection. 
J.  M.  Myers,  F.  H.  Madjid,  and  T.  T.  Wu,  Grant  N00014-89-J-1023  and 
Army  LABCOM  Contract  DAALQ2-89-K-0097;  Research  Unit  11. 

Antenna  arrays  are  potentially  useful  for  both  the  transmission  and  the  recep¬ 
tion  of  trains  of  electromagnetic  pulses.  For  example,  a  precisely  controlled  array  of 
antennas  driven  by  relativistic  klystron  amplifiers  (RKA’s)  can  potentially  launch 
trains  of  microwave  pulses  (say  N  pulses  per  train)  carrying  a  large  energy  flux — in 
Joules/m2 — to  strategic  distances.  Starting  from  information  given  to  it  about  a 
target,  an  effective  launcher  will:  1 )  deliver  adequate  cumulative  flux  over  N  pulses 
at  a  distance  as  great  as  the  target’s,  and  2)  aim  that  flux  to  illuminate  the  tar¬ 
get.  There  is  a  small  angle  by  which  the  microwave  beam  axis  can  miss  the  target 
and  still  deliver  adequate  flux.  This  angle  defines  an  allowable  aiming  tolerance. 
By  virtue  of  energy  conservation,  for  fixed  total  radiated  energy,  the  aiming  toler¬ 
ance  becomes  tighter  as  the  maximum  flux  is  increased.  This  trade-off  generates  a 
problem  in  aiming.  In  a  new  paper  [1],  this  aiming  problem  is  formulated  for  an 
RKA-driven  launcher  of  microwave  pulses,  and  an  approach  to  aiming  is  suggested 
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based  on  a  recent  study  of  quantum  measurements  [2]. 

For  the  RKA-launcher  (and  any  other  narrow-beam  device),  the  reward  for  the 
capacity  to  aim  the  beam  within  a  small  tolerance  6  is  to  deliver  a  flux  proportional 
to  9~2,  which  can  be  very  large.  The  difficulty  is  to  aim  the  beam  within  so  small  an 
angle.  In  testing  a  launcher,  it  is  desirable  to  aim  it  at  a  target  at  a  known  location. 
Knowledge  of  the  location,  however,  comes  with  an  uncertainty.  A  solution  has  been 
proposed  [1]  that  is  based  on  using  the  beam  itself  as  a  radar  to  reduce  uncertainty 
in  the  target  direction.  To  measure  the  target  direction  accurately  enough  to  reduce 
the  uncertainty,  one  can  augment  the  launcher  by  a  receiver  having  am  antenna  with 
aperture  diameter  comparable  to  that  of  the  launcher,  so  that  the  launcher  with  the 
aidjoining  receiver  acts  as  a  monopulse  radar  of  very  high  angular  precision.  The 
target  is  in  the  Fresnel  zone  and  not  the  far  zone  of  this  radar.  One  wamts  to  locate 
the  target  quickly  amd  accurately.  A  recent  study  of  trains  of  light  pulses  [2]  has 
shown  that  the  time  needed  to  make  a  quantum  measurement  of  given  accuracy  can 
be  reduced  by  varying  the  measurement  at  later  pulses  according  to  the  outcomes 
from  earlier  pulses.  This  suggests  that,  in  designing  a  radar  measurement  of  angle, 
one  look  for  a  parameter  to  vary  adaptively,  pulse  by  pulse,  guided  by  outcomes. 
The  appropriate  parameter  here  is  the  beamwidth,  to  be  varied  according  to  the 
dynamically  estimated  uncertainty  in  target  location. 
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IV.13  The  Launching  of  High-Power  Microwave  Pulses  from  a  Relativis¬ 
tic  Klystron  Amplifier.  T.  T.  Wu,  R.  W.  P.  King,  J.  M.  Myers,  and  H.- 
M.  Shen,  Grant  N00014-89-J-1023  and  Army  LABCOM  Contract  DAAL02- 
92-K-0036;  Research  Unit  11. 

A  new  project  initiated  in  1992  is  concerned  with  the  illumination  of  a  distant 
target  with  a  beam  of  electromagnetic  waves  that  is  characterized  by  both  high 
energy  and  high  frequency.  With  a  suitable  high-power  source  newly  available,  i.e., 
the  Relativistic  Klystron  Amplifier  (RKA),  the  problems  to  be  solved  relate  to  the 
development  and  design  of  an  efficient  system  for  transferring  power  from  the  source 
to  a  highly  directive,  accurately  steerable,  phased  array  of  radiating  elements.  The 
construction  of  an  apparatus  for  generating  and  projecting  a  very  narrow,  very- 
high-power  beam  of  electromagnetic  waves  that  is  both  adaptive  and  steerable  is 
not  a  simple  matter  of  design  using  familiar  components.  The  problems  involved 
require  am  innovative  approach  that  makes  use  of  a  wide  background  of  knowledge 
of  transmission  lines,  multi-mode  coaxial  and  single-mode  hollow  waveguides,  and 
antenna  theory  including  coupled  horns  and  parabolic  reflectors  as  elements  in  a 
steerable  phased  array. 

Answers  to  the  following  three  major  technical  questions  are  being  sought:  1) 
How  best  to  extract  the  EM  energy  from  an  RKA  and  transfer  it  efficiently  to  an 
antenna  system?  2)  How  to  design  the  antenna  system  so  that  it  will  transmit  the 
EM  energy  focused  into  a  narrow  beam?  3)  How  to  control  electrically  the  width 
and  direction  of  the  EM  beam  and  precisely  direct  it  to  a  target? 

Work  has  begun  on  the  first  question.  The  presently  suggested  way  to  trans¬ 
form  the  kinetic  energy  in  a  modulated  electron  beam  in  an  RKA  into  EM  energy 
is  by  means  of  a  coaxial- line  converter.  The  EM  field  generated  in  the  coaxial  line 
propagates  in  the  TEM  mode.  This  mode  is  not  suitable  for  transmitting  high 
power  to  the  elements  of  the  transmitting  antenna,  but  the  lowest  (H10)  modes  in 
rectangular  waveguides  are.  To  use  them,  it  is  necessary  to  have  a  suitable  matching 
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transformer  whose  purpose  is  to  transmit  the  TEM-mode  EM  field  from  the  coaxial 
line  into  the  H10-mode  of  a  number  of  rectangular  waveguides.  In  one  design  for  the 
transformer,  the  dimensions  of  the  coaxial  line  are  greatly  expanded  by  means  of  a 
biconical  section.  In  the  expanded  coaxial  line,  eight  fins  are  arranged  to  emerge 
from  the  inner  conductor  and  grow  gradually  until  they  reach  the  outer  conductor, 
and  so  divide  the  cross  section  of  the  coaxial  line  into  eight  sectors.  These  constitute 
separate  fan-shaped  waveguides  which  are  gradually  transformed  into  rectangular 
waveguides.  The  question  arises:  How  much  power  is  actually  transferred  to  the 
rectangular  waveguides?  Because  the  angular  distribution  of  the  electric  field  in  the 
TEM  mode  is  uniform,  the  presence  of  the  metal  fins  between  two  conductors  will 
cause  reflections.  To  find  out  how  much  energy  is  lost  in  the  reflection  and  conse¬ 
quent  mismatch,  it  is  necessary  to  calculate  the  composite  reflection  coefficient. 

The  analysis  of  the  RKA  coaxial-line  converter  has  been  initiated.  The  plan 
for  solving  this  problem  consists  of  two  steps:  solving  the  eigen-mode  problem  in  a 
large  coaxial  line  with  multi-fins  (CLMF)  of  constant  height,  and  then  varying  the 
height  of  the  fins  and  solving  the  coupling-mode  equation.  The  eigen-mode  in  a  large 
CLMF  with  fins  of  constant  height  has  been  cast  into  a  matrix  equation  which  has 
been  solved  numerically.  The  preliminary  results  show  that  the  eigenvalue  changes 
smoothly  with  the  height  of  the  fins.  When  the  height  is  small  or  zero,  the  eigenvalue 
is  also  small  or  zero,  like  the  eigenvalue  of  the  basic  TEM-mode  in  a  coaxial  line. 
When  the  fin  height  is  large  or  reaches  the  outer  conductor,  the  eigenvalue  is  close 
to  that  of  a  rectangular  waveguide.  The  analysis  of  the  RKA  converter  has  now 
advanced  to  the  second  step  of  formulating  the  coupling  equations  when  the  height 
of  the  fin  is  increasing  or  decreasing.  Of  primary  interest  is  the  evaluation  of  the 
reflected  mode  which  is  obtained  in  the  solution  of  the  coupling  equations  that 
determine  the  incident  and  reflected  modes. 
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IV. 14  Improved  Analysis  of  the  Resonant  Circular  Array.  G.  Fikioris, 
D.  K.  Freeman,  R.  W.  P.  King,  and  T.  T.  Wu,  Grant  N00014-89-J-1023, 
AF/RL  Contract  F19628-91-K-0020,  and  DOE  Grant  DE-FG02-84ER40158; 
Research  Unit  11. 

In  a  recent  paper  [1],  it  was  shown  theoretically  that  a  large  circular  array  of 
electrically  short  dipoles  may  possess  very  narrow  resonances.  The  array  consists  of 
N  identical,  perfectly  conducting,  parallel,  nonstaggered  elements;  only  element  1 
is  center-driven,  the  other  N  —  1  elements  axe  parasitic.  Subsequent  investigations 
[2]  have  shown  that  the  kernel  of  the  mth  phase- sequence  integral  equation  used 
in  [1]  does  not  take  into  account  the  inter-element  coupling  with  enough  accuracy 
for  a  resonant  circular  array.  The  formulation  in  [2]  of  the  mth  phase-sequence 
integral  equation  treats  self-  and  mutual  coupling  symmetrically,  by  averaging  the 
mutual  contribution  to  the  vector  potential  over  the  surface  of  the  dipole;  thus  the 
kernel  takes  the  form  of  a  double  integral.  The  imaginary  part  of  this  kernel — 
K\m\z)  (double  integral) — is,  under  suitable  conditions,  exponentially  small  in  IV, 
suggesting  the  existence  of  resonances  that  are  exponentially  narrow.  Furthermore, 
this  exponential  smallness  is  independent  of  fca;  i.e.,  when  K\m\z)  (double  integral) 
is  expanded  in  powers  of  fca,  each  term  is  exponentially  small. 

The  improved  kernel  developed  in  [2]  is  sufficiently  complicated  that  an  analyt¬ 
ical  study  of  the  two- term- theory  formulas  is  difficult  and  the  numerical  calculations 
involved  require  a  large  amount  of  computer  time.  More  recently,  an  alternative, 
simpler  version  of  the  kernel  in  [2]  has  been  proposed  [3]-[5]  which  possesses  many 
of  its  properties.  It  is  obtained  by  setting  a  =  0  in  the  self-term  of  the  imaginary 
part  of  the  kernel  used  in  (lj.  It  is  referred  to  as  the  “modified  kernel.”  The  leading 
term  in  the  asymptotic  expansion  of  Kjm\z)  (modified)  for  large  N  is  sufficiently 
simple  to  permit  an  analytical  study  of  the  two- term- theory  formulas.  The  conclu¬ 
sions  of  this  study  (supplemented  where  necessary  by  numerical  information)  have 
been  presented  as  a  sequence  of  twelve  basic  properties  relating  to  the  kernel,  the 
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conductances,  the  susceptances,  and  the  horizontal  and  vertical  field  patterns  [5]. 
In  addition,  two  refinements  to  the  modified  two- term  theory  are  described  that  are 
useful  for  numerical  calculations  whenever  high  accuracy  is  needed. 

The  analysis  above  is  concerned  with  lossless  elements.  For  the  case  of  highly 
conducting  elements — which  has  both  theoretical  and  practical  interest — a  new  for¬ 
mulation  of  the  theory  is  necessary.  In  the  mlh  phase- sequence  integral  equation, 
the  ohmic  losses  may  be  taken  into  account  with  a  small  term  proportional  to  the 
internal  resistance  per  unit  length  added  to  K\m\z).  By  modifying  Kjm\z)  in  this 
manner  throughout  the  formulation  of  the  two-term  theory,  a  new  solution  has  been 
obtained.  It  predicts  that  a  lossy  circular  array  as  seen  by  the  single  driven  voltage 
at  or  near  a  resonance  is  equivalent  to  a  series  RLC  circuit  in  which  R  =  Rtad  -f 
consists  of  two  resistances  in  series  (where  Rtkd  accounts  for  radiation  and  for 
the  ohmic  losses).  Both  of  these  are  limiting  factors  for  the  height  and  quality 
factor  of  the  resonance  curve.  The  situation  for  a  lossy  array  is  quite  different  from 
that  of  a  lossless  array:  If  Q(m,N)  is  the  quality  factor  in  the  lossless  case,  (a) 
Q  rises  rapidly  and  without  limit  as  the  number  N  of  elements  increases,  and  (b) 
Q  also  rises  rapidly  as  the  number  m  characterizing  the  phase-sequence  resonance 
increases.  In  the  case  of  a  lossy  array,  the  new  solution  predicts  that:  (a)  As  N 
increases,  Q  will  first  rise  rapidly  (while  ohmic  losses  are  negligible  and  radiation 
dominates)  and  then  decrease  linearly  (when  ohmic  losses  take  over  and  radiation 
becomes  small),  (b)  As  m  increases,  Q  will  first  rise  rapidly  (ohmic  losses  negligi¬ 
ble)  and  then  stay  constant  (ohmic  losses  dominant).  Thus,  the  new  theory  predicts 
that  it  is  no  longer  possible  to  obtain  resonances  that  are  as  narrow  as  one  wishes, 
as  in  the  lossless  case.  In  practical  arrays,  however,  only  moderately  narrow  reso¬ 
nances  are  required.  Actual  numerical  results  will  soon  be  obtained  using  this  new 
solution;  the  results  will  then  be  compared  with  the  experimental  data  (see  next 
topic)  to  determine  the  adequacy  of  the  theory. 
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In  order  to  verify  the  existence  of  the  theoretically  predicted,  extremely  sharp 
resonances  that  occur  in  a  circular  array  with  many  identical  elements  but  only 
one  driven,  experimental  investigations  of  such  an  array  have  been  carried  out.  Re¬ 
sults  from  the  first  experiment  have  been  reported  previously  [1];  they  showed  that 
resonance  does  occur.  They  also  demonstrated  how  sensitively  the  resonance  prop¬ 
erties  of  the  array  depend  on  the  precise  mechanical  implementation.  To  agree  with 
theory,  the  elements  should  be  identical,  perfectly  vertical,  and  positioned  exactly 
on  the  circle.  Any  slight  deviation  from  these  ideal  conditions  affects  the  mea¬ 
surements  significantly,  especially  the  resonant  admittance  values.  In  the  original 
experiment,  a  number  of  adjustments  to  the  array  were  needed  in  order  to  achieve 
close  agreement  between  the  measurements  and  the  end-corrected  theory  [2].  These 
adjustments  included  the  use  of  a  connecting  chain  to  improve  the  accuracy  of  the 
interelement  spacing,  the  addition  of  a  thin  disk  at  the  base  of  each  element  to  keep 
it  vertical,  and  low-resistance  conductive  cement  between  each  disk  and  the  ground 
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plane. 

To  improve  the  electrical  contact  between  the  individual  elements  and  the 
ground  plane,  and  to  remove  any  ambiguity  in  the  element  length  caused  by  the 
presence  of  the  thin  disk  at  the  base  of  each  element,  a  second  experiment  was 
designed  and  constructed  [3].  In  the  new  experiment,  the  choice  of  parameters  (i.e., 
the  height  h  and  radius  a  of  each  element  and  the  distance  d  between  adjacent 
elements)  is  based  on  the  results  of  the  modified  two-term  theory,  so  that  larger 
currents  around  the  array  should  be  obtained.  The  setup  consists  of  89  parasitic 
elements  and  one  driven  element,  equally  spaced  around  a  circle  of  40-in.  diame¬ 
ter.  The  single  driven  element  is  the  extension  of  the  inner  conductor  of  a  coaxial 
transmission  line.  The  parasitic  elements  are  thin-walled  brass  cylinders  that  are 
securely  fastened  (screwed)  onto  the  ground  plane  to  improve  the  electrical  contact 
with  the  ground  plane.  The  tolerances  in  the  location  of  the  holes  was  ±0.001  in. 
in  the  radial  positioning  and  ±0°0T'  in  the  angular  positioning.  The  monopoles 
were  machined  by  a  lathe  out  of  1/4- in.  brass  rods.  The  tolerance  in  the  heights  of 
the  parasitic  elements  was  ±0.001  in.  The  walls  were  made  quite  thin  so  that  the 
tubular  elements  of  the  theory  could  be  realized. 

The  purpose  of  the  experiment  is  to  calculate  the  driving-point  admittance 
of  the  array  from  swept-frequency  measurements  of  voltages  at  fixed  points  along 
the  transmission  line.  During  the  current  reporting  period,  different  measuring 
methods  have  been  developed  and  a  series  of  measurements  have  been  carried  out. 
Certain  difficulties  inherent  in  the  measurement  of  rapidly  varying  admittances 
were  encountered.  These  have  been  analyzed  and  several  approaches  to  overcoming 
these  difficulties  have  been  proposed.  Despite  the  difficulties,  it  has  been  shown 
that  narrow  resonances  occur  and  that  the  experimental  frequency-response  curves 
show  many  of  the  qualitative  features  predicted  by  theory.  Excellent  agreement 
has  been  obtained  between  the  experimentally  determined  resonant  frequencies  and 
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those  predicted  by  the  two-term  theory. 

The  finite  conductivity  of  the  elements  is  a  significant  limiting  factor  in  the 
large  values  of  current  that  can  be  attained  in  practice.  The  existing  theory  does  not 
take  the  finite  conductivity  into  account.  Investigations  are  in  progress  to  develop 
a  theory  that  includes  the  effect  of  ohmic  losses.  Once  this  has  been  accomplished, 
much  better  agreement  between  theory  and  experiment  is  expected  in  regard  to  the 
values  of  the  currents  at  resonance.  On  the  other  hand,  it  is  anticipated  that  ohmic 
losses  will  have  a  negligible  effect  in  shifting  the  theoretically  predicted  resonant 
frequencies. 
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IV.16  Properties  of  Closed-Loop  Arrays  of  Quantum  Mechanical  and 
Electromagnetic  Interactions.  D.  K.  Freeman  and  T.  T.  Wu,  Grant 
NG0014-89-J-1023  and  DOE  Grant  DE-FG02-84ER40158;  Research  Unit  11. 

This  research  project  has  now  been  completed.  The  doctoral  thesis  of  D. 
Kent  Freeman  [1]  has  been  defended  successfully  and  the  Ph.D.  degree  awarded 
to  Dr.  Freeman  in  June  1992.  The  purpose  of  this  study  has  been  to  trace  the 
phenomenon  of  extremely  narrow  resonances  in  closed- loop  arrays  from  its  simplest 
manifestation  in  arrays  of  Fermi  pseudopotentials  in  quantum  mechanics  to  its  cur¬ 
rent  status  as  a  means  of  designing  novel  electromagnetic  antenna  arrays.  The 
first  rigorous  investigation  of  resonances  in  closed-loop  arrays  was  carried  out  by 
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Grossmann  and  Wu  [2]  using  the  circular  array  of  uniformly  spaced  Fermi  pseudopo¬ 
tentials.  The  Fermi  pseudopotential  is  the  simplest  possible  nontrivial  scatterer  for 
the  Schrodinger  equation.  It  is  a  convenient  model  for  studying  array  phenomena 
that  occur  in  classical  electromagnetism  since  its  mathematical  formalism  is  consis¬ 
tent  and  explicitly  solvable.  This  simplification  is  justified  by  the  observation  that 
many  interesting  array  phenomena  depend  only  on  the  geometry  defined  by  the  po¬ 
sitions  of  the  scattering  centers  and  not  on  the  details  of  the  particular  interaction. 
It  was  found  that,  as  the  number  of  pseudopotentials  tends  to  infinity  while  the 
arc- length  spacing  is  held  constant,  such  an  array  admits  resonances  whose  widths 
tend  to  zero  exponentially.  This  result  was  the  foundation  for  the  present  work. 

Although  the  pseudopotential  provides  the  simplest  framework  for  studying 
the  resonance  phenomenon  of  closed  loops,  it  has  two  drawbacks.  First,  the  pseu¬ 
dopotential,  being  a  point  interaction,  lacks  structure;  thus,  it  is  impossible  to 
determine  what  effect  the  range  of  the  interaction  has  on  the  resonances.  To  ad¬ 
dress  this  question,  the  pseudopotential  has  been  generalized  to  an  interaction  with 
nonzero  range  by  studying  the  self-adjoint  integral-operator  potential  of  finite  rank. 
In  a  new  paper  [3] ,  it  is  shown  that  a  uniformly  spaced  infinite  linear  array  of  such 
potentials  can  admit  resonances  whose  widths  tend  to  zero  exponentially  as  the 
number  of  potentials  tends  to  infinity  while  the  arc- length  spacing  remains  con¬ 
stant.  Additionally,  it  is  demonstrated  that  the  pseudopotential  can  be  derived 
from  such  an  integral-operator  potential  as  am  appropriate  limit. 

The  second  drawback  of  the  pseudopotential  is  that  it  is  not  clear  if  results 
obtained  using  the  pseudopotential  are  applicable  to  antenna  arrays  in  classical 
electromagnetism.  Two  models  have  been  developed  for  studying  this  phenomenon 
directly  within  the  context  of  the  Maxwell  equations.  The  first  of  these  is  the 
pseudodipole,  which  was  introduced  by  Wu  [4]  as  the  simplest  possible  nontrivial 
scatterer  in  the  context  of  Maxwell’s  equations.  Like  the  pseudopotential,  the  pseu- 
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dodipole  is  a  point  interaction.  Therefore,  if  the  pseudodipole  is  to  approximate  a 
physical- dipole,  the  physical  dipole  should  be,  in  some  sense,  small.  Moreover,  in 
order  for  such  a  small  object  to  produce  a  significant  scattered  field,  one  expects 
that  it  must  be  resonant.  This  suggests  the  idea  of  talcing  a  center-loaded  cylindri¬ 
cal  dipole,  and  letting  its  size  tend  to  zero  while  simultaneously  adjusting  the  load 
reactance  in  order  to  keep  the  magnitude  of  the  scattered  field  constant.  In  a  recent 
paper  [5],  it  is  shown  that  for  an  incident  plane-wave,  this  procedure  results  in  a 
scattered  field  which  is  precisely  the  scattered  field  of  a  pseudodipole.  The  circular 
array  of  uniformly  spaced  pseudodipoles  has  also  been  studied  [6].  The  implica¬ 
tions  of  this  study  are  used  to  redefine  the  pseudodipole  in  a  more  mathematically 
consistent  manner,  as  well  as  to  further  elucidate  its  connection  with  the  physical 
cylindrical  dipole. 

The  second  of  the  two  electromagnetic  models  developed  in  [1]  is  a  new  formu¬ 
lation  for  studying  the  cylindrical  dipole  directly.  Perhaps  the  most  useful  method 
of  analysis  for  arrays  of  cylindrical  dipoles  is  the  “two- term  theory”  [7],  This  method 
has  the  advantage  of  producing  relatively  simple  closed-form  solutions  for  the  self- 
and  mutual  admittances  while  providing  remarkable  agreement  with  experimental 
results.  The  two-term  theory  is  not  accurate  enough,  however,  to  predict  correctly 
the  extremely  narrow  resonances  in  closed-loop  arrays  consisting  of  a  large  number 
of  dipoles.  This  problem  has  been  addressed  by  introducing  an  improved  kernel 
[8]  governing  the  interaction  between  cylindrical  dipoles.  This  kernel  is  inherently 
more  accurate  than  the  one  in  [7]  and  possesses  the  properties  required  to  predict 
correctly  the  narrow  resonances  in  closed-loop  arrays.  Finally,  in  order  to  gain  a 
better  understanding  of  why  the  two-term  theory  works  so  well,  a  new  formula¬ 
tion  of  this  theory  based  on  the  variational  principle  has  been  developed  [9].  It  is 
observed  that  the  two- term  theory  can  be  further  improved  by  incorporating  an 
end  correction  into  the  parasitic-dipole  current.  Combined  with  the  improved  ker- 
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nel,  this  new  formulation  provides  a  simple,  mathematically  sound  framework  for 
studying  resonances  in  closed-loop  arrays  of  cylindrical  dipoles. 
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Circular  closed-loop  arrays  of  parallel  dipoles  have  been  found  to  have  inter¬ 
esting  properties;  they  have  extremely  narrow  resonances  and  the  radiation  pattern 
is  also  extremely  narrow.  Unfortunately,  the  directivity  is  poor,  i.e.,  the  pattern  has 
many  equal  and  equally  spaced  lobes.  There  is  reason  to  think  that,  by  deforming 
these  arrays  from  the  circular  shape,  better  directivity  can  be  obtained.  The  search 
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for  closed-loop  arrays  with  high  directivity  is  a  challenging  task,  however,  since  there 
is  no  ready  theory  or  rule  to  follow.  It  is  extremely  difficult  to  develop  an  analytical 
method  for  determining  the  currents  on  the  elements  of  an  array  with  arbitrary 
shape  since  this  involves  solving  a  large  number  of  coupled  integral  equations. 

What  is  needed  is  a  qualitative  theory  that  can  determine  the  currents  ap¬ 
proximately  in  an  easy  but  physically  sound  way.  Based  on  an  understanding  of 
the  Yagi-type  parallel-dipole  array,  it  has  been  assumed  that  the  EM  wave  travels 
around  the  curved  line  guided  by  the  elements  with  their  currents  in  a  manner  char¬ 
acteristic  of  a  traveling  wave.  The  propagation  constant  is  assumed  to  be  the  same 
as  that  of  a  straight-line  array,  provided  the  radius  of  curvature  is  much  larger  than 
the  wavelength.  An  experiment  has  been  carried  out  to  verify  this  traveling-wave 
assumption  [1].  The  noncircular  array  is  obtained  by  displacing  only  a  small  num¬ 
ber  of  elements  from  the  original  circular  shape.  This  choice  of  shape  was  selected 
initially  simply  for  convenience  in  performing  the  experiment,  since  it  is  important 
to  keep  the  spacing  d  between  adjacent  elements  equal.  It  was  found  later  that 
such  a  shape  has  inherent  advantages  for  producing  higher  directivity.  The  detailed 
experimental  investigation  has  confirmed  that  a  traveling  wave  is  preserved  in  a 
curved  parallel-dipole  array.  The  small  energy  loss  or  amplitude  decay  associated 
with  a  curved  array  is  negligible  for  pattern  calculation. 

Once  the  currents  on  the  elements  have  been  determined,  calculating  the  ra¬ 
diation  pattern  is  straightforward.  However,  there  is  the  problem  of  how  to  search 
for  shapes  that  are  likely  to  produce  higher  directivity.  Since  there  are  innumerable 
possibilities,  it  is  hopeless  to  find  them  by  random  choice.  Consequently,  it  is  neces¬ 
sary  to  find  some  rules  which  can  guide  the  search.  The  traveling-wave  assumption 
described  above  has  been  used  to  study  a  variety  of  closed-loop  arrays,  and  guide¬ 
lines  concerning  radiation  from  a  traveling- wave  source  have  been  formulated  [1]. 
Based  on  these,  relations  have  been  established  between  individual  elements  (local 
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sources)  in  the  array  and  individual  lobes  (local  targets)  in  the  pattern.  It  is  found 
that  a  change  in  the  local  source  can  cause  a  change  in  its  related  local  target. 
Thus,  it  is  possible  to  locate  the  sources  of  the  individual  lobes  and  then  adjust 
these  sources  to  improve  the  pattern.  With  this  method,  patterns  have  been  im¬ 
proved  significantly;  the  main  lobe  has  been  increased  and  the  side  lobes  have  been 
substantially  reduced. 

Reference: 

1.  H.-M.  Shen,  “Experimental  Study  of  a  Noncircular  Array,”  SPIE  Proceedings, 

Vol.  1629,  Intense  Microwave  and  Particle  Beams  III,  281-297  (1992). 

IV.  18  Plasma  Waveguide:  A  Concept  to  Transfer  EM  Energy  in  Space. 
H.-M.  Shen,  Grant  N00014-89-J-1023  and  Army  LABCOM  Contract 
D A AL02-89- K-0097 ;  Research  Unit  11. 

It  has  been  shown  that  electromagnetic  (EM)  energy  originating  at  a  source 
in  upper  space  can  be  transferred  efficiently  through  a  plasma  waveguide  [1].  Such 
a  waveguide  consists  of  a  cylindrical  vacuum  core  surrounded  by  a  plasma  cladding 
medium,  which  can  be  generated  in  upper  space  by  a  hollow  laser  beam  [2].  An 
experiment  demonstrating  an  actual  95-meter  long  laser-generated  plasma  channel 
was  reported  in  1985  by  Livermore  National  Laboratory  [3].  The  analysis  and 
survey  show  that  the  plasma  waveguide  at  high  frequencies  (u>  »  up)  has  many 
excellent  performance  characteristics  for  transporting  a  microwave  pulse.  The  EM 
pulses  that  travel  in  such  a  waveguide  are  dispersion-free  and  have  a  speed  close  to 
that  of  light.  In  other  words,  the  EM  pulses  can  follow  the  laser-generated  plasma 
channel  and,  like  current  pulses  in  a  wire,  “move”  at  the  speed  of  light. 

In  the  first  model  of  the  plasma  waveguide  [1],  the  plasma  cladding  was  as¬ 
sumed  to  extend  transversely  to  infinity.  Therefore,  the  guided  modes  are  loss-free. 
If  the  cladding  is  finite,  there  will  be  an  energy  loss  from  the  outside  surface  of 
the  cladding.  To  study  this  effect,  an  analysis  of  the  plasma  waveguide  with  finite 
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thickness  was  carried  out  (4]  and  the  coefficient  of  the  energy  decay  calculated. 
It  was  found  that  most  of  the  properties  for  the  case  of  infinite  thickness  are  re¬ 
tained  for  the  plasma  waveguide  with  finite  thickness:  the  guided  modes  HEnm 
at  high  frequencies  are  split  into  two  frequency- independent  modes,  a  high-n  mode 
and  a  low-n  mode;  the  characteristic  equation  at  high  frequencies  is  also  split  into 
two  equations;  the  eigenvalue  x  (although  now  a  complex  number)  is  constant  and 
much  smaller  than  ka\  the  EM  pulse  propagates  with  constant  profile  and  shape 
and  at  the  speed  of  light;  and  the  loss  due  to  the  conductivity  of  the  plasma  and 
the  anisotropic  effect  due  to  the  earth’s  magnetic  field  are  small. 

The  main  concern  of  the  new  analysis  is  the  loss  from  the  outside  surface  of 
the  plasma  tube.  It  is  natural  to  suspect  the  loss  to  be  significant  since  the  field 
does  not  decrease  fast  enough  transversely.  For  instance,  the  field  on  the  surface 
at  r  =  2a  is  still  7%  of  the  maximum  field.  However,  the  results  of  this  analysis 
turn  out  to  be  quite  different.  For  the  case  of  finite  thickness,  the  eigenvalue  x  of 
the  characteristic  equation  becomes  a  complex  number.  The  imaginary  part  xt  is 
much  smaller  than  the  real  part  xR.  It  decreases  exponentially  with  the  thickness 
of  the  cladding.  Especially  for  high-frequency  operation — since  the  EM  wave  itself 
has  better  “directivity” — the  plasma  tube  has  the  ability  to  hold  most  of  the  energy 
within  the  tube,  i.e.,  the  loss  is  surprisingly  small. 

References: 

1.  H.-M.  Shen,  “Plasma  Waveguide:  A  Concept  to  Transfer  EM  Energy  in 
Space,”  J.  Appl.  Phys.  69(10),  6827-6835  (1991). 
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4.  H.-M.  Shen  and  H.-Y.  Pao,  “The  Plasma  Waveguide  with  a  Finite  Thickness 
of  Cladding,”  J.  Appl.  Phys.  70(11),  6653-6662  (1991). 
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ANNUAL  REPORT  OF 

PUBLICATIONS/PATENTS/PRESENTATIONS/HONORS 


a.  Papers  Submitted  to  Refereed  Journals  (and  not  yet  published) 

R.  W.  P.  King  and  S.  S.  Sandler,  “The  Detection  of  Dielectric  Spheres  Sub¬ 
merged  in  Water,”  IEEE  Trans.  Geosci.  &  Remote  Sensing  (accepted).  (Par¬ 
tial  support  from  the  Small  Business  Administration  under  the  SBIR  Program 
in  collaboration  with  Geo-Centers,  Inc.,  and  the  U.S.  Army  Medical  Branch.) 

R.  W.  P.  King,  “The  Propagation  of  a  Gaussian  Pulse  in  Sea  Water  and  Its 
Application  to  Remote  Sensing,”  IEEE  Trans.  Geosci.  &  Remote  Sensing. 
(Partial  support  from  Raytheon  Co.  under  Contract  N61533-90-C-0080  with 
ONR,  Annapolis  Detachment.) 

R.  W.  P.  King,  “The  Propagation  of  a  Low-Frequency  Burst  in  Sea  Water  and 
Its  Application  to  Remote  Sensing,”  IEEE  Trans.  Geosci.  &  Remote  Sensing. 
(Partial  support  from  Raytheon  Co.  under  Contract  N61533-90-C-0080  with 
ONR,  Annapolis  Detachment.) 

R.  W.  P.  King,  “Scattering  of  a  Cylinder  Submerged  in  the  Ocean:  The  Prop¬ 
agation  of  a  Low-Frequency  Rectangular  Pulse,”  Radio  Sci.  (Partial  support 
from  Raytheon  Co.  under  Contract  N61533-90-C-0080  with  ONR,  Annapolis 
Detachment.) 

R.  W.  P.  King  and  T.  T.  Wu,  “The  Propagation  of  a  Radar  Pulse  in  Sea 
Water,”  J.  Appl.  Pkys.  (Partial  support  from  Raytheon  Co.  under  Contract 
N61533-90-C-0080  with  ONR,  Annapolis  Detachment.) 

G.  Fikioris,  R.  W  .P.  King,  and  T.  T.  Wu,  “Improved  Analysis  of  the  Reso¬ 
nant  Circular  Array  of  Electrically  Short  Elements,”  J.  Appl.  Phys.  (Partial 
support  from  AF/RL  Contract  F19628-91-K-0020.) 

F.  H.  Madjid  and  J.  M.  Myers,  “Data-Engendered  Axioms  in  the  Quantum 
Theory  of  Optimal  Measurements,”  Annals  of  Physics.  (Partial  support  from 
Army  LABCOM  Contract  DAAL02-89-K-0097.) 

b.  Papers  Published  in  Refereed  Journals 

R.  W.  P.  King,  “Electromagnetic  Field  of  Dipoles  and  Patch  Antennas  on 
Microstrip,”  Radio.  Sci.  27(1),  71-78  (1992). 

R.  W.  P.  King,  “The  Circuit  Properties  and  Complete  Fields  of  Horizontal- 
Wire  Antennas  and  Arrays  Over  Earth  or  Sea,”  J.  Appl.  Phys.  71(3),  1499- 
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1508  (1992). 


M.  F.  Brown  and  J.  M.  Myers,  “Scattering  of  Electromagnetic  Missiles  by 
Convex  Metal  Objects,”  J.  Appl.  Phys.  70(11),  7176-7178  (1991).  (Partial 
support  from  Army  LABCOM  Contract  DAAL02-89-K-0097.) 

D.  K.  Freeman,  “The  Pseudodipole  as  the  Limit  of  a  Resonant  Cylindrical 
Dipole,”  J.  Math.  Phys.  32(7),  1961-1964  (1991).  (Partial  support  from  DOE 
Grant  DE-FG02-84ER40158.) 

D.  K.  Freeman,  “Generalization  of  the  Pseudopotential  to  an  Interaction  with 
Nonzero  Range,”  J.  Math.  Phys.  33(6),  2215-2227  (1992).  (Partial  support 
from  DOE  Grant  DE-FG02-84ER40158.) 

H.-M.  Shen  and  H.-Y.  Pao,  “The  Plasma  Waveguide  with  a  Finite  Thickness 
of  Cladding,”  J.  Appl.  Phys.  70(11),  6653-6662  (1991).  (Partied  support  from 
Army  LABCOM  Contract  DAAL02-89-K-0097.) 

d.  Books  (and  sections  thereof)  Published 

R.  W.  P.  King,  M.  Owens,  and  T.  T.  Wu,  Lateral  Electromagnetic  Waves: 
Theory  and  Applications  to  Communications,  Geophysical  Exploration,  and 
Remote  Sensing,  Springer- Verlag,  New  York  (1992). 

h.  Contributed  Presentations  at  Topical  or  Scientific/Technical  Society 
Conferences 

J.  M.  Myers  and  H.-M.  Shen,  “Experiments  Using  Laser- TViggered  Electrical 
Discharges  to  Generate  Electromagnetic  Missiles,”  SPIE  Proceedings,  Vol. 
1629,  Intense  Microwave  and  Particle  Beams  III,  274-280  (1992).  (Partial 
support  from  Army  LABCOM  Contract  DAAL02-89-K-0097.) 

H.-M.  Shen,  “Experimental  Study  of  a  Noncircular  Array,”  SPIE  Proceedings, 
Vol.  1629,  Intense  Microwave  and  Particle  Beams  III,  281-297  (1992).  (Par¬ 
tial  support  from  AF/RL  Contract  F19628-91-K-0020  and  Army  LABCOM 
Contract  DAAL02-89-K-0097.) 

F.  H.  Madjid,  J.  M.  Myers,  and  T.  T.  Wu,  “Lessons  for  Pulsed-Array  Radar 
from  Quantum  Light  Detection,”  SPIE  Proceedings,  Vol.  1629,  Intense  Mi¬ 
crowave  and  Particle  Beams  III,  368-376  (1992).  (Partial  support  from  Army 
LABCOM  Contract  DAAL02-89-K-0097.) 

i.  Honors/Awards/Prizes 

R.  W.  P.  King,  1991  Distinguished  Service  Award,  from  the  IEEE  Antennas 
and  Propagation  Society. 


102 
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for  the  Year  Ending  30  June  1992 

Mr.  V.  Houdzoumis. 


103 


104 


V.  SIGNIFICANT  ACCOMPLISHMENTS  REPORT 


V.l  Competitive  Analog  Neural  Networks  for  Translation-Invariant  Fea¬ 
ture  Detection.  F.R.  Waugh  and  R.M.  Westervelt;  Research  Unit  5,  SOLID 
STATE  ELECTRONICS. 

Competition  among  a  group  of  neurons  is  common  in  pattern-recognizing  neu¬ 
ral  networks.  Competitive  winner-take-all  networks  can  mimic  associative  memory 
behavior  if  each  neuron  performs  a  dot-product  of  the  network  input  with  the  pat¬ 
tern  that  the  neuron  represents.  Much  less  common  are  networks  of  groups  of 
competing  neurons,  which  are  functionally  similar  to  networks  of  interconnected 
associative  memories.  Fukushima’s  neocognitron  is  a  feed-forward  example  of  such 
a  network. 

We  have  studied  the  dynamics  of  recurrent  networks  of  interconnected  groups 
of  competing  analog  neurons.  We  have  shown  that,  with  symmetric  interconnections 
and  discrete-time,  parallel  updating,  these  networks  are  stable  when  the  slopes  of  the 
neuron  input-output  transfer  functions  are  sufficiently  small.  When  configured  in 
a  two-dimensional,  locally- connected  architecture  (see  accompanying  figure),  these 
networks  can  perform  translation-invariant  feature  detection  on  visual  images.  Each 
neuron  in  a  competing  group  represents  a  fragment  of  the  feature  to  be  detected,  and 
neurons  receive  inputs  both  from  a  localized  region  of  a  pixel  array  and  from  neurons 
in  neighboring  clusters.  Interconnections  between  clusters  enable  the  networks  to 
reconstruct  complete  features  from  their  constituent  fragments,  with  translation 
invariance  on  the  fragment  length  scale. 

Competition  can  be  implemented  naturally  in  analog  electronics  by  forcing 
competing  neurons  to  share  a  limited  resource,  such  as  the  current  through  a  com¬ 
mon  wire.  Moreover,  competition  can  greatly  reduce  the  density  of  neural  intercon¬ 
nections  in  visual  processing  applications,  because  neurons  communicate  informa¬ 
tion  about  whole  features  rather  than  individual  pixels. 
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COMPETITIVE  ANALOG  NEURAL  NETWORKS  FOR 
TRANSLATION-INVARIANT  FEATURE  DETECTION 


architecture 
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array 


2D  competitive 
network 


(a)  competing  neurons  represent 
pattern  fragments 


(b)  connections  between  neurons  store 
relations  among  pattern  fragments 
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V.2  A  Novel  Resonant  Array  of  Parallel  Dipoles  with  Unique  and  Useful 

Properties.  R.  W.  P.  King,  T.  T.  Wu,  D.  K.  Freeman,  G.  Fikioris,  and 
H.-M.  Shen;  Research  Unit  11,  ELECTROMAGNETIC  PHENOMENA. 

A  theoretical  and  experimental  study  has  been  made  of  the  remarkable  prop¬ 
erties  of  a  circular  array  of  a  large  number  N  of  parallel  dipoles  with  only  one 
element  driven.  The  length  and  radius  of  each  tubular  element  and  the  distance 
between  adjacent  elements  have  been  determined  analytically  for  the  particular  sets 
of  values  for  which  the  entire  array  resonates.  A  large  number  of  possible  modes  ex¬ 
ists.  For  each,  the  currents  in  the  elements  attain  extremely  large  amplitudes  with 
phase  relations  determined  by  the  particular  mode.  The  driving-point  impedance  of 
the  single  driven  element  is  a  pure  resistance  at  each  resonance.  The  experimental 
array  consists  of  90  elements  and  the  predictions  of  the  theory  have  been  verified 
over  a  range  of  frequencies  that  includes  a  sequence  of  many  resonant  modes. 

The  field  pattern  of  the  highest  mode  (in  which  the  phases  of  the  currents 
in  adjacent  elements  differ  by  exactly  180°)  consists  of  a  circle  of  90  pencil- like 
spikes,  each  with  extremely  narrow  beam  width  in  both  vertical  and  horizontal 
planes.  This  mode  has  an  interesting  application  as  a  beacon.  When  rotated  at 
a  desired  frequency — e.g.,  1  rpm — the  array  sends  out  90  pulses  per  minute  in  all 
directions.  If  a  unidirectional  pattern  is  desired,  the  90  elements  can  be  rearranged 
in  a  serpentine  manner  with  the  condition  of  resonance  preserved.  An  optimum 
shape  and  design  leading  to  a  true  supergain  pattern  remains  to  be  determined.  In 
another  important  modification,  a  single  omnidirectional  pancake-like  pattern  with 
an  extremely  narrow  vertical  beam  width  can  be  obtained.  When  located  over  the 
earth  or  sea,  this  generates  a  powerful  omnidirectional  surface  wave  with  virtually 
no  sky  wave.  Application  to  a  ground-wave  over-the-horizon  radar  is  indicated. 


107 


Significant  Accomplishments  for  FY1992 

Research  Unit  11  on  ELECTROMAGNETIC  PHENOMENA 
Harvard  University — JSEP  Grant  N00014-89-J-1023 
Principal  Investigator:  Professor  Tai  Tsun  Wu 

A  Novel  Resonant  Array  of  Parallel  Dipoles 
with  Unique  and  Useful  Properties 


Highly  accurate  formulas  have  been  developed  and  verified  ex¬ 
perimentally  for  the  properties  of  a  closed  loop  of  a  large  num¬ 
ber  of  dipoles  of  which  only  on  or  two  are  driven. 


Resonant 

array 


Unusual  Properties  in  Resonant  Mode: 

••  Driving-point  impedance:  Very  small  pure  resistance. 

••  Field  patterns  governed  by  shape  of  loop  and  excitation. 

•  N  radial  spikes;  rotate  for  use  as  beacon. 

•  Single,  highly  directional  beam. 

•  Omnidirectional,  very  flat,  pancake-like  pattern;  de¬ 
ploy  array  in  space  for  communication;  erect  over 
earth  or  sea  for  ground-wave  communication,  or  over- 
the-horizon  radar  in  (^J.000  km  range. 


